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ABSTRACT 

Recent  gas  flow  studies  on  linear  systems  noted  the 
importance  of  the  interaction  of  molecular  streaming  and 
inertial  effects.  This  led  to  speculation  regarding  the 
nature  of  these  effects  in  radial  systems.  As  a  result,  it 
was  decided  to  examine  radial  flow  both  experimentally  and 
theoretically . 

Radial  gas  flow  theory  has  been  modified  to  account  for 
viscous,  inertial,  and  molecular  streaming  effects.  In  order 
to  carry  out  the  experimental  verification  of  the  proposed 
theory,  flow  tests  were  conducted  on  three  Indiana  limestone 
samples  measuring  approximately  one  foot  in  diameter  and  one 
to  one  and  a  half  inches  in  thickness.  A  radial  flow  cell 
was  specially  designed,'  built,  and  tested  for  the  purpose. 

The  resulting  data,  when  examined  both  graphically  and 
numerically,  indicated  that  all  three  effects  are  important, 
and  that  the  suggested  modification  of  the  graphical  procedure 
leads  to  more  realistic  parameter  estimations.  It  has  also 
been  found  that  the  Klinkenberg  graphical  method  is  inapplic¬ 
able  under  certain  conditions,  and  that  the  manner  in  which  a 
test  sample  is  confined  during  a  laboratory  flow  test  is  a 
factor  of  considerable  importance. 
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INTRODUCTION 


The  phenomenon  of  gas  flow  through  porous  media  has 
been  a  subject  of  investigation  for  quite  some  time.  Numerous 
investigators  have  carried  out  experiments  to  verify  Darcy’s 
equation  for  viscous  flow  and  For chheimer ’ s  modification 
of  it  for  high  velocity  flow.  The  phenomenon  of  gas-slipp¬ 
age,*  studied  by  Klinkenberg,  has  also  been  experimentally 
verified  in  many  cases. 

Considerable  research  has  been  done  on  these  subjects 
at  the  University  of  Alberta  and  many  interesting  results 
obtained.  However,  the  majority  of  this  work  and  that  done 
elsewhere  has  been  on  linear  systems. 

In  general,  the  flow  in  natural  gas  reservoirs  is 
roughly  radial.  Hence,  the  question  arises  whether  the 
results  obtained  on  linear  systems  may  be  applied  to  the 
actual  radial  reservoirs  by  a  mere  transformation  of  geo¬ 
metry,  or  whether  some  factors  hitherto  unaccounted  for  come 
into  play. 

Because  of  the  relative  scarcity  of  results  on  radial 
systems,  this  study  was  undertaken.  From  a  theoretical 
standpoint  this  necessitated  a  reconsideration  of  the  basic 
theory  of  the  flow  of  gases  in  porous  media,  especially  in 
light  of  the  recent  studies  on  linear  systems  which  empha- 


*  the  terms  gas-slippage  and  molecular  streaming  are  used 
interchangeably  in  this  work. 
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size  the  importance  of  molecular  streaming.  For 
out  the  experimental  verifications,  new  equipment 
and  assembled  to  carry  out  the  gas  flow  tests. 


carrying 
was  designed 
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LITERATURE  REVIEW 

# 


The  available  literature  on  the  subject  of  fluid  flow 
through  porous  media  is  quite  abundant.  To  encompass  the 
various  aspects  of  this  topic,  the  present  section  is  sub¬ 
classified  to  review  the  work  done  in  the  areas  of  viscous 
flow,  vis  co-inert ial,  flow  and  their  experimental  verifications. 

Viscous  Flow 

The  experiments  of  Henry  Darcy  (1856),  concerning  the 
flow  of  water  through  sand  filters,  was  the  first  scientific 
study  of  fluid  flow  through  porous  media.  As  a  result  of 
this  and  further  investigations  the  following  relationship, 
known  as  Darcy's  law,  was  proposed  to  described  steady 
horizontal  flow  of  fluids  in  the  viscous  region. 


IE  =  Hi 

dx  k 


(1) 


M.  King  Hubbert  1  ^  has  shown  the  theoretical  soundness 
of  this  experimental  proposition  by  deriving  it  from  the 
fundamental  Navier-S t okes  equation  of  motion  for  viscous 
fluids . 

Although  equation  (1)  holds  in  case  of  liquid  flow,  it 

(2) 

was  shown  by  Klinkenbergv  '  that  the  permeability  constant 
as  determined  with  gase,s  is  dependent  upon  the  nature  of 
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. 


4 


the  gas,  and  is  approximately  a  linear  function  of  the 
reciprocal  mean  pressure.  This  phenomenon  was  attributed 
to  molecular  streaming.  In  view  of  this,  the  permeability 
k  in  equation  (1)  should  be  replaced  by  an  apparent  perme¬ 
ability  k^  defined  as 


k  =  k(l  +  -) 
a  p 


(2) 


Under  the  assumptions  of  isothermal,  viscous  flow  of 

(3) 

gas  in  a  homogeneous  medium,  Al-Hussainy  e_t  a_l.  studied 

the  effect  of  variations  of  pressure-dependent  gas  viscosity 
and  compressibility.  For  steady  radial  gas  flow  in  the 
viscous  region  they  proposed  the  relation 


TTkh  T  [m  (p  )  -  m  (p  )  ] 

sc  r  e  w 


s  c 


T  p  In  — 
sc  r 


w 


(3) 


where 


m(p) 


P 

y(p)  Z(p) 


(4) 


dp 
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Vi s co - Ine r t i al  Flow 

It  was  observed  that,  when  the  fluid  flow  rate  was 
continuously  increased,  Darcy's  law  did  not  adequately  des¬ 
cribe  the  flow  behaviour  beyond  a  certain  flow  rate.  To 

(4) 

describe  this  high  velocity  flow,  Forchheimer  suggested  a 
modification  of  Darcy's  law  and  proposed  an  equation 

=  aq  +  bq2  (5) 


where  a  and  b  are  thought  to  be  constants. 

This  model  was  obtained  on  the  basis  of  empirical  reasoning 
and  curve  fitting  techniques  applied  to  various  experimental 
data  when  available. 

Adopting  an  approach  similar  to  that  of  M,  King 
Hubbert, Chwyl^^  derived  Forchheimer ' s  equation  from  the 
fundamental  Navier-Stokes  equation,  thereby  putting  it  on  a 
sound  theoretical  basis. 

Green  and  Duwez^^  demonstrated  the  validity  of  the 
equation 

“  dx  =  aPq  +  FBPq2 

for  describing  high  velocity  linear  flow.  On  the  basis 
of  this,  they  defined  the  Reynold's  number  for  porous  medium 


as 


6 


Vq 

ay 


(7) 


where 


a 


1 

k 


and  F  =  the  inertial  resistance  coefficient . 

D 

Using  the  modified  gas  law  as  an  equation  of  state, 

/  Q  \ 

Greenberg  and  Wegerv  integrated  equation  (6)  to  yield 


Mp2)Mg  _  +  F  (G) 

2L  (  ZRT)  yG  B^y' 


(8) 


where  G  =  pq  =  gas  mass  velocity 

2 

and  g  =  32.17  lbm- f t / lb f -s e c  . 

Defining  Reynold's  number  as 


Nr 


Dvp 


yx 


n 


and  the  Fanning  friction  factor  as 


dp  _  f v  p 
dr  =  2gDXn 


and  utilizing  the  friction-factor  vs.  Reynolds  number  charts 
for  Brownell  and  Katz,  the  following  integrated  forms  for 
radial  systems  for  laminar  and  turbulent  regions  were  given 


udiun  bn&  ■xenlmsl  ioi  em938Y3  iBibfli 
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by 


Elenbaas  and  Katz 


.(9) 


(P 


Pw} 


=  0.000812 


y  ZTW 


2  ,  n-m 
D  <p  hG 


In 


w 


(9) 


(P 


Pw} 


=  0.00000202 


W2TZf 

D  nh2G 


w 


(10) 


where  W  =  mass  flow  rate  in  Ibs/sec 

D  =  average  diameter  of  particles  in  ft 
(j)  =  porosity  of  the  bed 
h  =  bed  thickness  in  ft 
G  =  gass  gravity 

and  n  and  m  are  the  exponents  to  be  determined  as  suggested 
in  their  paper . 

Combining  the  modified  gas  law  with  the  radial  equivalent 
of  equation  (6),  Tek  e t  a_l .  presented  the  following  inte¬ 

grated  equation, 


P  “  P 
e  w 


1424  yZT  Q  In ( — — ) 
0  r 

_ w_ 

hk 


+ 


3.1602  x  10 


-12 


F  CTZ 
B 


(id 


where  y  =  viscosity  in  centipoise  evaluated  at  arithmetic 
mean  temperature  and  pressure 

Z  =  compressibility  factor  evaluated  at  arithmetic 
mean  temperature  and  pressure 


. 


i  inss  nt  ^3 jgooalv  *  14 
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T  =  average  temperature  of  the  flowing  gas,  °R 
Qq  =  volumetric  flow  rate,  MSCF/day 
k  =  permeability  in  millidarcies 
and  Fg  =  inertial  resistance  coefficient,  ft  ^  . 

This  so  far  is  the  most  widely  accepted  model  for  steady 
radial  gas  flow. 

Assuming  porous  media  to  be  made  up  of  bundles  of 
capillary  tubes,  with  orifice  plates  spaced  throughout  the 
tubes  at  a  distance  equal  to  the  tube  diameter,  Blick^^ 
proposed  a  capillary  orifice  model.  Applying  the  force 
momentum  balance  and  modified  gas  law,  he  derived  the  follow¬ 
ing  equation  for  describing  steady  linear  isothermal  gas  flow, 


2  2  2 

(pl  -  p 2 ) M6  4> 

2LG2ZRT 


(12) 


where 


and 


16 

Nr 

e 


for  viscous  flow 


4.0  log  0  (Nr ^  /cj)  -  0.4,  for  turbulent  regime 


CD  =  drag  coefficient  of  orifice  plate 
G  =  pv 

6  =  average  pore  diameter,  ft. 


and 


porosity  fraction. 
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As  a  result  of  this  development,  he  was  able  to  develop  the 
following  theoretical  expression  for  the  inertial  resistance 
coefficient  : 


PC 


D 


B 


2A4>26 


(13) 


where  A  -  1.0. 

(12) 

Extending  Blick'sv  technique  to  radial  systems, 


C  raf  t  on 
flow  : 


(13) 


suggested  the  following  equation  for  radial  gas 


(P 


2 

e 


18546 


Q  p  H  ZT 

w  w 


khG 


+  5.41  x  10 


2  2 

-10  FBQWPWZT 
h2G 


w 


3  r 


(14) 


where  =  well  bore  flow  rate,  MCF/day 

3 

p^  =  well  bore  fluid  density,  lbm/ft  . 

Experimental  Investigations 

Many  experimental  investigations  have  been  carried  out 
to  verify  the  validity  of  Darcy’s  law  for  gas  flow  through 
porous  media  and  Klinkenberg ' s  slippage  phenomenon.  Notable 


s  dcoM  . nonsmonariq  9*s„  .3  »' J f  =•  dii^i  IA  bus  *lb*K  ,'uoioq 
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among  these  are  the  works  of  Estes  and  Fulton,  Rose, 

Grunberg  and  Nissan,  and  Calhoun  and  Yuster . 

On  the  basis  of  flow  experiments  conducted  on  porous 

(18) 

sintered  metals,  Greenberg  and  Weger ^  J  reported  that  per¬ 
meability,  while  not  a  function  of  pressure,  varied  inversely 
with  temperature.  They  concluded  that  the  inertial  resistance 
coefficient  remained  constant  for  pressures  up  to  2000  psia 

in  the  temperature  range  15°C  to  60°C. 

( 19 ) 

Stewart  and  Owens  studied  the  extent  to  which 

slippage  and  inertial  effects  could  exist  simultaneously,  by 
performing  tests  on  heterogeneous  limestones.  Though  at 
very  large  Reynold's  numbers  the  slippage  effects  disappeared, 
it  was  found  that  both  slippage  and  inertial  effects  were 


present  at  high  flow  rates. 

Dranchuk  and  Sadiq^®) 


have  shown  the  need  for  proper 


delineation  of  the  viscous  and  visco-iner tial  regions,  so 
that  the  Klinkenberg  plot  of  equations  (2)  and  the  visco- 
iner tial  plot  of  equation  (8)  give  consistent  permeability 


values  . 

(21) 

Hamilton  considered  the  slippage  to  be  significant 

at  low  mean  pressures.  Using  equation  (8)  he  obtained  the 
values  of  k  and  from  the  straight  line  portion  of  the  visco 

inertial  plot.  The  downwarping  of  the  plot  at  lower  flow 
rates  was  considered  to  be  due  to  slippage  effects. 

Investigating  the  magnitude  of  resulting  errors,  when 
average  values  for  viscosity  and  compressibility  of  the  fluid 


■ 
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(22) 

are  used,  Mackett  concluded  that  the  values  of  the 

apparent  gas  permeability  calculated  using  these  properties 
did  not  differ  appreciably  from  the  values  obtained  using 
integrals  to  express  the  isothermal  variation  of  these  pro¬ 
perties  . 

While  re-examining  the  data  of  previous  investigators, 

(2  3) 

Kolada  observed  that  invariably  the  permeability  obtained 

from  a  Klinkenberg  plot  was  lower  than  that  given  by  a  standard 
vis co-iner t ial  plot.  This  discrepancy  was  attributed  to 
molecular  streaming  effects  which  had  been  neglected  in  the 
vis co-iner tial  flow  theory.  Considering  molecular  streaming 
for  flow  under  both  low  and  high  pressure  gradients,  he  deve¬ 
loped  a  modified  form  of  the  vis co-iner tial  gas  flow  theory 
for  linear  systems  and  also  presented  its  experimental  verifi¬ 
cation. 

(24) 

Based  on  the  above  work,  Dranchuk  and  Kolada  have 

suggested  a  modified  technique  for  interpreting  laboratory 
data  obtained  in  gas  flow  tests  conducted  on  linear  porous 


media  . 
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STEADY  RADIAL  GAS  FLOW  WITH  SLIPPAGE 


(25) 

The  work  of  Dranchuk  and  Kolada  on  linear  systems 

has  established  the  need  for  incorporating  molecular  stream¬ 
ing  effects  into  the  basic  vis  co-inert ial  gas  flow  theory. 

As  seen  in  the  linear  case,  the  molecular  streaming 
effects  being  important,  their  consideration  in  the  radial 
vis  co-inert ial  gas  flow  theory  is  warranted.  Examination 

/  Q  Z'  \  /  Q  T  \ 

of  the  available  models  of  Tek  e_t  a_l.  and  of  Crafton^ 

shows  that  they  do  not  allow  for  this  effect.  However,  it 

appears  that  they  could  be  modified  in  order  to  account  for 

(28) 

gas  slippage.  In  fact,  the  Tek  £_t  al .  model  was  so  modi 

(  29  ) 

fied.  Modification  of  Crafton's  model  though  planned, 


was  not  completed  due  to  shortage  of  time. 


A  Modified  Radial  Gas  Flow  Equation 

The  modification  of  the  Tek  e_t  a  1 .  model  may  be 

explained  by  means  of  the  following  reasoning. 

Steady  radial  horizontal  gas  flow  under  the  conditions 
of  no  gas  slippage  is  described  by  the  relationship 


_  dp  =  yq 

dr  k  B 


pq !  q 


(15) 


In  this  equation,  the  absolute  value  of  ’q*  is  employed  be- 


13 


cause  'q?  is  conventionally  taken  as  positive  when  flow  is 
in  the  direction  of  increasing  'rf  and  negative  when  it  is 
in  the  opposite  direction. 

Using  the  modified  gas  law  together  with  average  values 
of  viscosity,  compressibility  and  temperature,  the  above 
equation  on  integration  yields  the  following  form  which  is 

(31) 

essentially  the  same  as  that  suggested  by  Tek  et_  al . 

except  that  the  sign  convention  for  'q'  or  'Qq’  ^as  been  used. 


2  2 

<’e  -  *v>  ‘  -  [' 


1424  yZTQ0  ln(^) 


w 


hk 


+ 


3.1602  x  10  12  F  GTZQ  |Q  |(^-  - 

w  e 


(16) 


Under  the  conditions  of  gas  slippage,  equation  (15) 
has  to  be  modified  so  that, 


1e  =  Ha 

dr  k 

a 


+ 


Fgpq  ! q 


Multiplying  through  by  p  yields 


d£  =  ]£££ 

p  dr  k 

a 


V2q 


(17) 


(18) 


rf:  ,  9i  l  ij  i  sqm  jJ  bnfi  73il  i  c  .2  ■  ri  q  a  *j  o  d  ,i(Jiso3aiv  3to 
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By  virtue  of  the  modified  gas  law, 


pM 

ZRT 


(19) 


and 


P 


0 


P0M 


RT 


0 


(20) 


Furthermore,  for  plane  radial  flow, 


q 


Q  =  Q 

A  2  7Trh 


(21) 


an  d 


p0Q0  =  PQ 


(22) 


Using  the  average  values  of  viscosity  and  compressibi¬ 
lity  evaluated  at  the  arithmetic  mean  temperature  and  pres¬ 
sure  and  substituting  equations  (19),  (20),  (21),  and  (22)  in¬ 
to  equation  (18)  yields, 


M  p  dp 
ZRT  dr 


P0Q0MU 

RT^k  (  2  TTrh ) 
0  a 


P0M  2 

+  F„(~-) 


B  'RT 


0 


V  V 

(  2  TTrh )  2 


(23) 


ol  ,  noun  3 Ji J  turi 


' 

.  &  (81)  no  ups  oU 
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o  r 


uzTPqQo 

k  T  ( 2  irrh) 
a  0 


+  F 


(P0M) 2  ZRT  Q0  I  % 


B  'RT 


0 


M 


(  2irrh ) 


(24) 


but 


k  =  k(l  +  -) 
a  p 


(25) 


Hence 


pzTPqQq 

k  ( 1  +  -)  TA(2TTrh) 
P  0 


+  F 


PqM  2 

(^T~) 


B  RT 


0 


ZRT  Q0  I  Q0  I 
M  (  2  TTrh  )  2 


(26) 


Defining 


C 


1 


P?.Tp0Q0 

kT  (2uh) 


(27) 


and 


q  - 


(28) 


C 


2 


P0M  2 


B  RT 


0 


ZRT  %  I Q0 
M 

(  2  TTh ) 


2 


equation  (26)  can  be  concisely  written  as 


(1  +  ~)  r 


which  may  also  be  written  as 


C  r  +  C  (1  +  ^) 

-1 _ £ _ E_ 

(1  +  j)r2 


(29) 


(30) 


Equation  (30)  is  a  non-linear  ordinary  differential 
equation  for  which  there  is  as  yet  no  known  solution.  How- 
ever,  if  —  is  sufficiently  small,  this  equation  may  be 
approximated  by  means  of  the  equation 


Clr  +  C  2 


(31) 


which  is  separable.  Equation  (31)  may  be  written  as 


-  p(l  +  ~)  d  p  =  (-—  +  — |)dr 

p  r 


(32) 


as  n  j.  ihw  yioaionoo  sd  ns  (d£)  nol:fsup9 
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and  integrated  between  the  limits, 


r  =  r  ,  p  =  p 
w  w 


r  =  r  , 
e 


p  =  p 


to  give 


(P2  -  P2)  r 

-  I— y  “  +  b(Pe  '  p„)]  =  C1 

w 


+  c2(7-  -  -T- >  (33) 

w  e 


When  multiplied  throughout  by  -2.0  this  equation  becomes 


2  2  re 

(P  -  p  )  +  2b(p  -  p  )  =  -  2  [  C  ln( — ) 
vte  *w  re  rw  1  r 

w 


+  C2(^-f)]  (34) 

w  e 


If  the  customary  field  units  of 


p  =  centipoise,  h,  r  ,  rg  =  ft, 

k  =  millidarcy,  Qq  =  MSCF/day, 


T 


0 


5  2  0  0  R 


14.7  p s ia 
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and  G  =  gas  gravity,  where  gravity  of  air  =  1.0  are 
employed,  this  equation  takes  the  form, 


(p2  “  P2)  +  2b(p  -  p  )  =  -  [ 

e  w  e  w 


1424  pZTQ0  In 
_ w 

hk 


(35) 


3.1602  x  10  12  FbGTZQ0  |Qq  |(^-  - 

+  - 3 - 2 - S_] 


It  is  interesting  to  note  that,  although  equation  (35) 
is  but  an  approximate  solution  to  the  differential  equation 
describing  flow  where  slippage,  viscous, and  inertial  effects 
are  important,  neglecting  inertial  effects  causes  the  equation 
to  collapse  to  the  form' 


r 

1424  yZTQ0  ln(^) 

(p2  -  p2)  +  2b(pe  -  Pw)  -  -  [ - hk - - 1  (36) 


which  is  in  fact  the  rigorous  solution  for  that  case. 
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USE  OF  THE  PROPOSED  EQUATIONS  FOR 
PARAMETER  ESTIMATIONS 


The  forms  of  both  the  equations  (35)  and  (36)  are  such 
that  they  are  quite  amenable  to  graphical  treatment.  How¬ 
ever,  in  case  of  equation  (35),  one  of  the  three  parameters 
must  be  either  known  or  assumed  in  order  to  solve  for  the 
other  two.  For  example,  if  T  b ’  is  known  a  priori,  'k'  and 

'F  '  may  be  evaluated. 

B 

In  case  of  flow  towards  the  well  bore,  rearrangement  of 

2  2 

equation  (36)  after  multiplication  by  k/(p  -  p  )  gives 

0  w 


k(l  +  3) 
P 


r 

1424  uZTQq  ln(^) 

w 


)h 


(37) 


It  can  be  recognized  that  the  right  hand  side  of  the  above 
equation  is  in  fact  the  apparent  permeability  k  .  Hence,  it 

Si 

along  with  —  may  be  used  to  construct  the  Klinkenberg  plot 
P 

for  the  radial  flow  case,  from  which  'k'  and  'b'  may  be  deter¬ 
mined  . 

(  32 ) 

As  observed  by  Dranchuk  and  Sadiq,  the  Klinkenberg 

plot  gives  meaningful  results  provided  the  data  used  are  in 
the  viscous  region.  This  presents  the  difficulty  of  de¬ 
lineating  the  viscous  and  vis  co-inert ial  flow  regions.  By 
employing  the  Reynold's  number  criterion,  this  difficulty  can 


ffiria  ad  ago  3l 
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be  overcome.  Since  the  Reynold’s  number  for  a  porous  medium 
is  defined  as  the  ratio  of  inertial  effects  to  viscous  effects, 
it  can  be  obtained  by  taking  the  ratio  of  the  second  term 
to  the  first  term  on  the  right  hand  side  of  equation  (35). 

It  may  be  observed  that  equation  (35)  can  be  rearranged 
to  yield 


p^)+2b(p  -p)] 

w  e  *w 

- -  - j 

1424  yZTQ0 


[ 


1«£> 

w 


(38) 


+ 


2.2192  x  10~15  FBG(— 

B  r 

w 


1  l^o- 

— )  ( - — ) 

r  y 
e 


When  the  value  of  1  b ’  obtained  in  the  previous  plot  is 
used,  equation  (38)  suggests  a  straight  line  plot  whose  inter¬ 
cept  and  slope  should  yield  the  values  of  ’k’  and  ’F  '  res- 

D 

pectively.  On  the  other  hand,  equation  (16)  may  be  placed 
in  the  following  form 


2  2 

[(pe  ~  Pw)h] 


[-^  + 


2.2192  x  10  15  F  G (— 

B  r 


w 


Iq 


0 


(39) 


1424  yZTQ 


0 


y 


which  may  also  be  plotted  in  order  to  estimate  'k'  and  'F  '. 
However,  in  this  case,  the  line  passing  through  high  flow  rate 
data  points  which  presumably  lie  on  a  straight  line  must  be 
extrapolated . 


' 
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In  view  of  the  fact  that  the  plot  made  employing  equa¬ 
tion  (38)  is  the  radial  flow  equivalent  of  Kolada  's 
modified  vis  co-inert ial  plot,  and  the  plot  made  employing 
equation  (39)  corresponds  to  his  conventional  vis  co-inert ial 
plot;  these  two  plots  are  referred  to  as  the  modified  and 
conventional  vis  co-inert ial  plots  respectively. 

Since  these  graphical  methods  are  very  similar  to  those 

(34) 

employed  by  Kolada  for  the  linear  case,  they  pose  similar 

problems.  For  example,  the  curvature  in  the  conventional 

vis  co-inert ial  plot  makes  the  delineation  of  the  straight  line 

portion  of  the  curve  difficult  if  not  impossible.  On  the  other 

hand,  determination  of  ’b’  from  the  Klinkenberg  plot  and  'k' 

and  'F  '  from  the  modified  vis  co-inert ial  plot  requires  data- 
B 

splitting.  This  in  turn  requires  not  only  a  large  number  of 

data  points  but  approximately  the  same  number  of  points  for 

each  plot  if  the  ’k’  values  obtained  from  the  two  plots  are  to 

be  consistent.  Furthermore,  since  data  splitting  is  based  on 

the  Reynold’s  number  criterion,  the  delineation  of  the  viscous 

region  is  not  only  a  trial  and  error  procedure,  but  it  is  also 

dependent  upon  the  Reynold's  number  one  chooses. 

From  this  it  is  obvious  that  a  numerical  technique, 

which  utilizes  all  of  the  data  simultaneously  and  does  not 

require  the  delineation  of  flow  regimes,  would  be  preferable 

to  the  graphical  methods.  One  such  technique  is  that  of 

(  35  ) 

quasilinearization  suggested  by  Kolada. 
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EXPERIMENTAL  APPARATUS 


In  order  that  plane  radial  gas  flow  through  porous  media 
be  exmained  experimentally,  samples  of  Indiana  limestone, 
resembling  large  grind  stones,  were  obtained  from  Imperial 
Oil  Production  Research  Laboratory  in  Calgary.  These  samples 
were  about  one  foot  in  diameter  and  between  one  and  two 
inches  thick. 

Prior  to  flow  tests,  the  porosity  of  these  samples  was 
determined  using  a  specially  constructed  Boyle's  law  porosi- 
meter.  This  essentially  consisted  of  a  core  holder  and  an 
expansion  chamber  connected  to  gas  supply  and  a  vacuum  pump 
as  schematically  shown  in  Figure  1.  By  calibrating  with 
solid  aluminium  blanks  of  known  dimensions,  the  gross  grain 
volume  of  the  cores  was  .estimated  and  porosity  determined. 

In  order  to  conduct  flow  tests  on  the  available  samples, 
a  special  core  holder( as  described  below,  was  designed  and 
cons  true  ted . 

Core  Holder 

A  mild  steel  cylinder  having  an  outer  diameter  of  18 
inches,  a  height  of  3.862  inches,  and  a  wall  thickness  of 
2.50  inches  was  placed  between  two  mild  steel  blind  flanges 
which  were  22  inches  in  diameter  and  2.75  inches  thick  as 
shown  in  Figure  2.  A  seal  between  the  three  pieces  was 
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FIGURE  1  .  FLOW  DIAGRAM  OF  THE  EQUIPMENT  FOR  POROSITY  DETERMINATION 


RADIAL  FLOW  CELL 
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effected  by  means  of  ’O'  ring  grooves  and  'O'  rings  positioned 
in  the  top  and  bottom  of  the  cylinder.  A  close  fitting  pis¬ 
ton  was  placed  inside  the  cylinder.  The  piston  was  constructed 
of  mild  steel,  was  1.778  inches  thick,  and  had  an  ’O’  ring 
and  two  back  up  rings  in  its  periphery.  This  piston  could 

be  raised  or  lowered  hydraulically,  by  means  of  an  Owattona 

3 

series  y21-l  model  !C'  hydraulic  pump.  Four  g-  inch  holes  were 
drilled  around  the  circumference  of  the  cylinder  to  serve  as 
gas  inlets,  while  a  central  hole  in  the  top  flange  was  pro¬ 
vided  to  serve  as  the  gas  outlet.  In  preparation  for  a  run, 
the  sample  was  placed  on  the  piston  between  two  ^  inch  thick 
natural  rubber  gaskets  cut  to  the  exact  shape  of  the  rock 
samples.  The  blind  flanges  were  then  bolted  together  by  means 
of  ten  equiangular ly  placed  one  inch  thick  bolts.  The  piston 
was  raised  until  the  rock  sample  was  confined  between  it  and 
the  upper  flange. 

Apparatus  Assembly 

The  apparatus  assembly  is  schematically  shown  in  Figure 
3.  The  four  gas  inlets  were  connected  to  a  common  manifold 
constructed  of  one  inch  pipes  .  The  manifold  in  turn  was 
connected  to  the  gas  supply  cylinder  via  regulators  and 

which  were  used  for  controlling  low  and  high  upstream 
pressures  respectively.  The  gas  outlet  was  connected  to  a 
soap  film  burette  and  wet  test  meter  which  were  used  for  mea¬ 
suring  low  and  high  flow  rates  respectively.  The  wet  test 
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FIGURE  3.  FLOW  DIAGRAM  FOR  ASSEMBLED  EQUIPMENT  FOR  GAS  FLOW 


meter  was  in  turn  connected  to  an  ’Angus'  automatic  strip 
chart  recorder  with  an  event  marker.  While  the  mercury  mano¬ 
meter  MM^  and  Heise  pressure  gauge  PG^  recorded  the  upstream 
pressure,  the  mercury  manometer  MM^  and  the  pressure  gauge 
PG^  recorded  the  downstream  pressure.  Copper-cons tantan 
thermocouples  TC-^  and  TC^  were  used  to  measure  the  inlet  and 
outlet  temperatures  of  the  flowing  gas.  Pressure  gauge  PG^ 
indicated  the  confinement  pressure  on  the  sample. 
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ROCK  SAMPLES  AND  FLOW  TEST  METHODS 


A  total  of  three  samples  of  Indiana  limestone  were 
studied.  Two  of  these  were  obtained  from  the  Imperial  Oil 
Research  Laboratory  in  Calgary.  The  third  was  cut  in  the 
workshop  of  the  Department  of  Chemical  and  Petroleum  Engineer¬ 
ing,  University  of  Alberta,  from  Indiana  limestone  slabs 
supplied  by  Sharp  Bros.  Cut  Stone  Company,  Hamilton,  Ontario. 

The  samples  obtained  from  the  Imperial  Oil  Laboratory 
were  somewhat  more  heterogeneous  than  the  one  prepared  at  the 
University,  since  they  had  been  acid  treated  in  the  vicinity 
of  the  inner  radius  before  being  subjected  to  the  flow  studies 
described  herein. 

Employing  a  combination  of  one  of  two  different  methods 
of  maintaining  confining  pressure  on  the  sample  and  one  of 
two  different  ways  of  changing  the  flow  rate,  a  total  of 
four  different  types  of  flow  tests  were  conducted.  These 
tests  were  conducted  in  the  following  manner. 

(1)  The  confining  pressure  was  held  constant  throughout 
the  test  and  flow  rate  was  changed  by  raising  the  up¬ 
stream  pressure  while  holding  the  downstream  pressure 
at  atmospheric . 

(2)  The  sample  was  confined  as  in  (1)  above  and  the  flow 
rate  was  increased  in  stepwise  fashion  by  decreasing 
the  downstream  pressure  while  the  upstream  pressure  was 
held  constant  at  a  certain  fixed  value. 
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(3)  The  flow  rate  variation  technique  was  identical  to  that 
in  (1)  above  but  the  confining  pressure  was  increased 
with  each  flow  rate  increase  so  that  the  difference  be¬ 
tween  the  confining  pressure  and  the  arithmetic  mean  flow¬ 
ing  pressure  was  held  constant  at  some  preselected  value. 

A  combination  of  the  flow  rate  variation  technique  employed 
in  (2)  and  the  confining  pressure  technique  employed  in 
(3)  was  used. 


(4) 


. 
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EXPERIMENTAL  PROCEDURE 


The  pressure  gauges,  wet  test  meter,  and  thermocouples 
were  calibrated  before  connecting,  them  to  the  experimental 
apparatus.  All  the  fittings  and  connections  were  pressure 
tested  to  ensure  a  leak-proof  assembly. 

The  outer  diameter,  inner  diameter,  and  thickness  were 
obtained  for  each  sample  by  taking  the  average  of  ten  caliper 
measurements.  Bulk  volumes  were  calculated  using  these  aver¬ 
age  values.  Porosity  was  determined  using  the  equipment 
previously  described. 

The  sample  was  properly  mounted  inside  the  sample  holder 
and  subjected  to  the  desired  confining  pressure. 

Gas  from  the  supply  cylinder  was  flowed  through  the 
appropriate  regulator  at  the  desired  pressure,  and  was  allowed 
to  flow  for  a  sufficient  time  to  ensure  stabilized  conditions. 
Inlet  and  outlet  pressures  and  temperatures  and  corresponding 
flow  rate  measurements  were  taken  using  the  appropriate 
devices.  Barometric  pressure,  temperature  of  the  flowing 
gas  prior  to  entering  the  flow  devices,  and  the  pressure  shown 
by  water  manometer  ( ' WM ’  in  Figure  3)  were  also  recorded. 
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TREATMENT  OF  DATA 


Volumetric  flow  rates  of  Nitrogen  at  the  recorded 
temperatures  and  pressures  (given  by  manometer  f WM ’ )  were 
converted  to  the  reference  conditions  of  one  international 
atmosphere  (14.696  Psia)  and  520°R,  using  the  modified  gas 
law.  Compressibility  factors  at  required  temperatures  were 
determined  by  five-point  Lagrangian  interpolation  of  the 

/  O  £  \ 

curve  fitted  data  of  Hilsenrath  e_t  al . 

Viscosity  values  at  the  arithmetic  mean  temperature 

and  pressure  of  the  flowing  gas  were  obtained  using  the  follow- 

(37) 

ing  relationship  suggested  by  Kestin  and  Wang. 


y  =  1778  x  10_7(1  +  8.958  x  10~4(p-l) 

+  6.120  x  10"7(p-l)2  +  3.997  x  10_8(p-l)3) 

+  4.55  x  10~  7 (T-25 )  (40) 

where  y  is  in  poise 

p  is  in  international  atmospheres 
and  Tisin°C. 

A  computer  program  was  written  to  evaluate  the  expres¬ 
sions  necessary  for  making  the  Klinkenberg  permeability  plot. 
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the  conventional  vis co-iner t ial  plot,  and  the  modified  visco- 
inertial  plot.  Listings  of  this  program  are  given  in  Appendix 
C.  A  program  was  also  written  to  perform  a  linear  least 
squares  fit  of  the  data  wherever  needed. 

The  parameters  k,  b,  and  F  were  evaluated  by  applying 

D 

q ua sil inea r i z a t ion  to  equation  (35).  The  technique  of  qua¬ 
silinearization  has  been  discussed  in  detail  by  Kolada. 

The  salient  points  of  quasilinearization  as  applied  in  the 
present  case,  are  given  in  Appendix  A.  A  sample  of  the  com¬ 
puter  program,  written  to  effect  the  quasilinearization  solu¬ 
tion,  is  given  in  Appendix  C. 

Quasilinearization  being  an  iterative  process,  starting 
from  given  initial  guesses,  it  is  important  to  have  good 
initial  guesses  to  assure  convergence. 

Permeability  was  approximated  by  that  obtained  from 
the  lowest  flow  rate  observation.  This  value  of  permeability 
together  with  the  following  correlation  suggested  by  Heid 
et  al.^9)  was  used  for  estimating  the  value  of  f  b  T 


b 


0 . 777 


-0.39 

k 


(41) 


where  b  is  in  international  atmospheres  and  k  is  in  milli- 

darcies .  Approximation  for  the  inertial  resistance  coefficient 

J  v  -  ,  (40) 

was  obtained  using  the  correlation  suggested  by  Kolada. 
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F 


B 


7.56  x  108 

^1.67  kl . 12 


(42) 


where  k  =  millidarcies  and  (j)  =  fractional  porosity. 

In  summary,  the  following  procedure  was  adopted  for 
parameter  estimations  by  various  methods. 

(1)  Values  of  k,b,  and  were  estimated  using  qua¬ 
silinearization. 

(2)  On  the  basis  of  these  values  the  highest  flow  rate, 
for  which  Reynold’s  number  was  less  than  0.01, 

was  determined. 

(3)  The  data  points,  for  which  the  flow  rates  were 
less  than  that  determined  in  (2)  above,  were 
plotted  on  a  Klinkenberg  plot  and  values  of  k  and 
b  determined  using  a  linear  least  squares  fit. 

(4)  Using  the  value  of  'b'  determined  in  (3)  above, 

the  modified  vis  co-inert ial  plot  was  made  to  yield 
the  values  of  k  and  F  from  a  linear  least  squares 
fit.  -  B 

(5)  For  comparison,  the  conventional  vis  co-inert ial 
plot  was  also  made. 
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RESULTS 


The  physical  properties  of  the  samples  studied  were 
as  shown  in  Table  1. 

Sample  111  was  the  first  to  be  studied.  A  Klinkenberg 
plot  for  it  is  presented  in  Figure  4.  Departure  from  linearity, 
marking  the  transition  to  vis  co-inert ial  flow,  appears  to 
start  at  flow  rates  higher  than  0.0937  MSCF/day.  This  checks 
with  the  Reynold's  number  criterion  of  Nr^  =  0.01.  Using  the 
value  of  'b '  obtained  from  this  plot,  the  modified  visco- 
inertial  plot  was  made  as  shown  in  Figure  5.  A  conventional 
vis  co-inert ial  plot  also  appears  in  the  same  figure. 

Starting  with  the  initial  guesses  of 

k  =  4.679  md 

b  =  6.255  psia 

and 

Fg  =  1.679  x  108  ft"1 

q uas i 1 ine a r i z a t i on  was  applied  using  equation  (35)  and  a 

-4 

convergence  criterion  of  £  =  10  .  Finally  the  parameters 

obtained  by  means  of  quasilinearization  were  used  to  generate 
data  which  were  plotted  in  the  form  of  a  conventional  visco- 
inertial  plot  as  shown  in  Figure  5.  A  sample  calculation  for 
this  method  of  data  generation  is  given  in  Appendix  B. 


PHYSICAL  PROPERTIES  OF  THE  SAMPLE 
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FIGURE  5.  VISCO-INERTIAL  FLOW  PLOTS  SAMPLE  1 
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The  parameter  estimations  obtained  graphically  and 
numerically  are  summarized  in  Table  2. 

TABLE  2 

A  SUMMARY  OF  RESULTS  FOR  SAMPLE  111 


Me  t ho  d 

k 

(md) 

b 

(psia) 

fb. 

(ft 

-1) 

Klinkenberg  Plot 

2  .  704 

12.495 

Conventional  Vis co-Iner t ial 

^  9 

Plot 

3 . 102 

8.198 

x  10 

Modified  Vis co-Inertial 

q 

Plot 

2.727 

12 . 49  5* 

6.287 

x  lO9 

Quasilinearization 

2 .719 

12.600 

6.142 

x  109 

*  in  the  above  table  and  succeeding  tables,  an  underlined 

value  of  b  indicates  that  the  value  was  not  determined  from 
that  plot  but  in  fact  used  in  its  construction. 


Adopting  the  same  procedure  as  in  the  case  of  sample  111, 
of  varying  the  flow  rates  by  holding  the  downstream  pressure 
at  atmospheric  and  increasing  the  upstream  pressure  stepwise, 
run  1  on  sample  112  was  recorded.  A  Klinkenberg  plot  and 
vis  co_inertial  plots  for  this  run  were  constructed  as  shown 
in  Figures  6  and  7  respectively.  From  Figure  6  it  appears 
that  data  points  corresponding  to  flow  rates  higher  than 
0.0209  MS CF / day  fall  outside  the  viscous  region. 
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FIGURE  6.  KLINKENBERG  PLOT  SAMPLE  112  RUN  1 
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FIGURE  7.  VISCO-INERTIAL  FLOW  PLOTS  SAMPLE  112  RUN  1 
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Using 


k  -  16.617  md 

b  =  3.815  psia 

and 

Fg  =  9  .235  x  108  ft"1 

as  initial  guesses,  five  iterations  were  required  for  the 

quasilinearization  to  attain  the  convergence  criterion  of 

-4 

e  =  10  . 

The  results  obtained  by  these  various  methods  are 
tabulated  in  Table  3. 

Runs  2  and  3  for  the  same  sample  were  taken  at  constant 
upstream  pressures  of  100  psig  and  80  psig  respectively 
with  flow  rates  being  gradually  increased  by  decreasing  the 
downstream  pressure.  As  can  be  seen  from  Figure  8,  this 
procedure  resulted  in  peculiar  Klinkenberg  plots.  The  nature 
of  these  plots  is  such  that  they  cannot  be  subjected  to  con¬ 
ventional  interpretation  in  so  far  as  the  prediction  of 
permeability  and  slippage  coefficient  is  concerned. 

Though  at  first  glance  these  plots  may  appear  to  be 
suspect,  the  following  reasoning  suggests  that  this  is  in 
fact  what  should  be  expected.  When  the  downstream  pressure 
is  constant  and  upstream  pressure  is  successively  increased, 

decreases  while  the  flow  rate  increases,  resulting  in  lower 
P 
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FIGURE  8.  KL1NKENBERG  PLOT  SAMPLE  112  RUNS  2  AND  3 


TABLE  3 


A  SUMMARY  OF  RESULTS 

FOR  RUN  1 

SAMPLE 

112 

Me  tho  d 

k 

(md) 

b 

(psia) 

fb 

(ft 

X) 

Klinkenberg  Plot 

12.34 

5.00 

Conventional  Vis co-Inert ial 
Plot 

13.78 

5.94  x 

109 

Modified  Vis co-Inert ial 

Plot 

12 .26 

5.00 

5.62  x 

109 

Quasilinearization  of 
Equation  (35) 

12.31 

4.57 

5  .  55  x 

109 

apparent  permeabilities.  However,  for  a  constant  upstream 

pressure  when  the  downstream  pressure  is  decreased,  _ 

P 

increases.  But,  since  the  flow  rate  also  increases,  it  should 
result  in  lower  apparent  permeabilities.  This  being  precisely 
what  is  seen  in  Figure  8,  the  given  plots  are  quite  in  order. 

Figure  (9)  shows  the  Klinkenberg  plots  for  the  three 
runs  on  sample  112.  As  stated  earlier,  runs  2  and  3  were 
recorded  by  imposing  back  pressure.  The  difficulty  experi¬ 
enced  when  subjecting  these  curves  to  the  usual  interpreta- 
tion  technique  suggests  that  Sadiq's(41)  recommendation 
concerning  the  use  of  back  pressure  to  bring  the  points  up  to 
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P  Ipsia 

FIGURE  9.  KLINKENBERG  PLOTS  SAMPLE  112  RUNS  1,  2  AND  3 


the  straight  line  may  not  be  valid. 

Despite  the  fact  that  these  runs  produced  anomolous 
Klinkenberg  plots,  the  resulting  vi s c o— iner t ia 1  plots  appear 
to  be  normal  in  every  respect,  as  may  be  seen  from  Figures 
10  and  11.  It  should  be  noted  that,  as  a  result  of  anomolous 
Klinkenberg  plots,  parameter  estimation  by  graphical  means 
is  not  possible  in  these  cases  and  hence  the  modified  visco- 
inertial  plots  cannot  be  constructed.  However,  the  modified 
vis co -iner t ial  plots  shown  in  these  two  figures  were  con¬ 
structed  using  the  value  of  'b'  obtained  from  quasilineariza¬ 
tion  and  are  presented  merely  to  show  that  the  resulting 
behaviour  was  normal. 

A  comparison  of  results  obtained  for  runs  1,  2,  and  3 
on  sample  112  is  presented  in  Table  4. 

Run  1  on  sample  2,  recorded  under  constant  downstream 
pressure,  showed  a  very  unusual  behaviour  when  plotted  on  a 
Klinkenberg  plot  (Figure  12).  Contrary  to  expectation  the 
higher  flow  rate  data  points  showed  an  increase  in  apparent 
permeability  with  increasing  mean  pressure.  Though  the  tenta¬ 
tive  viscous  data  gave  a  permeability  of  0.639  md  and  a 
slippage  coefficient  of  7.158  psia,  neither  the  vis  co-inert ial 
plots  (Figure  13)  nor  quas i 1 ine ar i z a t i on  gave  meaningful 
values.  In  fact,  both  of  them  gave  negative  inertial  resistance 
coefficient  values  which  have  no  known  physical  meaning. 


Run  2  on  the  same  sample,  which  was  recorded  with  a 
constant  upstream  pressure  of  200  psig  exhibited  behaviour 
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FIGURE  10.  VISCO-INERTIAL  FLOW  PLOTS  SAMPLE  112  RUN  2 
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RESULTS  OF  VARIOUS  RUNS  ON  SAMPLE  112 
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FIGURE  12.  KLINKENBERG  PLOTS  SAMPLE  2  RUNS  1  AND  2 
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FIGURE  13.  ViSCO-iNERTIAL  FLOW  PLOTS  SAMPLE  2  RUN  1 


51 


similar  to  that  observed  in  runs  2  and  3  on  sample  112,  as 
can  be  seen  from  Figure  12. 

Although  the  conventional  vis co-inertial  plot  for  this 
run  was  smooth  its  slope  was  such  that  extrapolation  was 
impossible.  On  the  other  hand,  the  modified  vis co-inertial 
plot  yielded  the  expected  straight  line  as  is  evident  from 
Figure  14 . 

The  application  of  quasilinearization  to  the  data  for 
this  run  resulted  in  the  following  estimates. 


k  =  0.73  md 

b  =  35.87  psia 

and 

FB  =  2  .  71  x  1011  ft"1 

It  should  be  noted  that  because  of  the  nature  of  the 
Klinkenberg  plot,  parameter  estimation  by  graphical  means 
alone  was  not  possible  and  that  the  modified  plot  in  Figure  1-f 
was  constructed  using  the  value  for  'b  obtained  numerically. 
Furthermore,  the  modified  plots  indicated  that 


k  =  0 . 74  md 


and 
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FIGURE  14.  ViSCO-INERTIAL  FLOW  PLOTS  SAMPLE  2  RUN 
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which  merely  shows  that  behaviour  was  normal. 

The  abnormal  behaviour  in  Run  1  was  thought  to  be  an 
effect  of  confining  pressure.  It  was  reasoned  that  increases 
in  mean  flowing  pressures  resulted  in  decreases  in  net  confin¬ 
ing  pressures,*  which  in  turn  resulted  in  increases  in  apparent 
permeability.  The  fact  that  this  effect  was  not  very  apparent 
in  Run  2  was  thought  to  be  due  to  the  difference  in  test 
p  ro  cedure  . 

In  order  that  this  interpretation  of  the  behaviour  might 
be  evaluated,  run  3  was  made.  In  this  case  the  net  confining 
pressure  was  held  constant  at  400  psig  and  the  downstream 
pressure  was  held  at  atmospheric.  To  further  elucidate  this 
behaviour,  run  4  was  conducted  by  maintaining  a  constant  net 
confining  pressure  of  400  psig  and  holding  the  upstream  pres¬ 
sure  constant  at  200  psig. 

As  may  be  seen  from  Figure  15,  the  Klinkenberg  plots 
for  both  these  runs  exhibit  normal  behaviour.  Figure  15A 
shows  the  results  of  run  4  and  part  of  run  3  on  an  enlarged 
scale.  Vis  co-inert ial  plots  for  these  runs  are  shown  in 
Figures  16  and  17  again  indicating  normal  behaviour.  From 
these  plots  it  was  concluded  that  the  reasoning  with  regards 
to  variations  in  net  confining  pressure  was  essentially  correct. 

The  results  of  these  various  parameter  estimation  tech¬ 
niques  as  applied  to  runs  3  and  4  are  summarized  in  Table  5. 

*  Net  confining  pressure  =  Imposed  confining  pressure 
-  arithmetic  mean  flowing  pressure. 
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FIGURE  16.  VISCO-INERTIAL  FLOW  PLOTS  SAMPLE  2  RUN  3 
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FIGURE  17.  VISCO-INERTIAL  FLOW  PLOTS  SAMPLE  2  RUN  4 


The  observed  basic  flow  data  and  quasilinearization 
results  for  all  the  runs  on  the  various  samples  are  tabulated 
in  Appendix  D.  Comparison  of  the  observed  function  values 
and  those  calculated  using  the  parameters  obtained  from  qua¬ 
silinearization  together  with  the  relative  percentage  error 
between  the  two  also  appear  at  the  end  of  each  table  in  the 


same  Appendix. 


SUMMARY  OF  RESULTS  FOR  RUNS  3  AND  4  ON  SAMPLE 
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DISCUSSION  OF  RESULTS 


The  modified  form  of  the  vis  co-inert ial  flow  equation  is 
quite  amenable  to  graphical  treatment  provided  the  slippage 
coefficient,  b,  is  known  a  priori. 

Permeability  values  obtained  from  the  Klinkenberg  plot 
check  very  well  with  those  predicted  from  the  modified  visco- 
inertial  plot  and  quasilinearization  as  can  be  seen  from 
Tables  2,  4  and  5.  Also  the  slippage  coefficient  values 
derived  from  the  Klinkenberg  plot  and  quasilinearization  show 
good  agreement. 

Graphical  treatments  require  a  large  number  of  observa¬ 
tions  since  the  data  have  to  be  split  to  adequately  cover 
each  region.  Moreover,  as  stated  earlier,  they  have  another 
drawback  in  that  the  data  splitting  has  to  be  done  on  the 
basis  of  a  Reynold's  number,  the  calculations  of  which  requires 
the  knowledge  of  k  and  F  .  Since  these  parameters  happen  to 
be  the  ones  being  estimated,  a  trial  and  error  procedure  is 
required . 

On  the  other  hand,  the  numerical  technique  of  qua¬ 
silinearization,  since  it  utilizes  all  of  the  data  simul¬ 
taneously  and  is  rapid,  is  preferable  if  access  to  a  digital 
computer  is  available. 

This  method  does  have  one  serious  shortcoming  in  that 
it  uses  the  data  irrespective  of  their  quality.  As  a  result 
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poor  or  what  may  be  called  "wild"  data  points  cannot  be 
detected  in  this  method,  as  readily  as  in  the  graphical 
approach  as  may  be  seen  from  Figure  11. 

Permeability  values  obtained  from  conventional  visco- 
inertial  plots,  in  general,  do  not  correspond  to  those 
obtained  from  Klinkenberg  extrapolations.  Furthermore,  when 
the  slippage  coefficient  is  high  it  is  difficult,  if  not 
impossible,  to  extrapolate  the  assumed  unaffected  data. 

This  is  well  illustrated  in  Figures  5,  10  and  14. 

Even  when  the  slippage  coefficient  is  small,  the  con¬ 
ventional  plot  does  not  yield  inertial  resistance  coefficient 
values  corresponding  to  those  derived  from  quasilinearization, 
although  permeability  values  may  show  good  agreement.  In 
this  respect,  the  better  agreement  in  results  with  other 
techniques  shows  the  superiority  of  the  modified  plotting 
technique . 

One  of  the  most  important  observations  made  during  this 
study  was  that  the  shape  of  the  Klinkenberg  plot  is  not 
unique  but  very  much  depends  on  the  manner  in  which  a  flow 
experiment  is  conducted.  Figures  8,  9,  12  and  15  bring 

this  out  quite  clearly  . 

As  shown  by  the  results  of  Run  1  on  Sample  2  (Figure  12) 
uncontrolled  confining  pressures  could  lead  to  unrealistic 
data  and  consequent  serious  misinterpretations.  Hence,  the 
control  of  confining  pressure,  at  least  in  certain  cases, 
could  very  well  be  an  important  factor  of  consideration. 
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In  all  of  the  runs,  taken  while  holding  the  upstream 
pressure  constant,  the  observed  points  were  outside  the 
viscous  region.  This  is  to  be  expected  since, at  higher 
pressures,  even  a  small  pressure  drop  results  in  a  large 
flow  rate  . 

Permeability  values  predicted  from  the  flow  experiments 
conducted  at  constant  upstream  pressure  compare  well  with 
those  obtained  from  the  experiment  with  constant  downstream 
pressure.  However,  the  slippage  and  inertial  resistance 
coefficients  show  different  values  in  each  case.  It  is  not 
readily  clear  whether  this  is  due  to  some  shortcomings  in 
either  the  experimental  or  interpretation  technique  or  whether 
it  is  a  manifestation  of  some  hitherto  unknown  effect. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

As  a  result  of  this  study,  the  following  conclusions 

were  drawn: 

(1)  The  approximate  solution,  to  the  suggested  differential 
equation,  appears  to  describe  reasonably  well,  the  radial 
gas  flow  through  porous  media. 

(2)  The  conventional  vis co-inertial  plot  does  not  seem  to 
be  adequate  for  estimation  of  permeability  and  inertial 
resistance  coefficient  values. 

(3)  The  modified  vis  co-inert ial  plot  and  the  Klinkenberg 
plot  may  be  used  to  estimate  values  of  permeability, 
slippage  and  inertial  resistance  coefficients.  However, 
the  accuracy  of  these  appears  to  be  somewhat  lower  than 
that  realized  using  quasilinearization. 

(4)  Gas  slippage  seems  to  affect  the  flow  behaviour  even  at 
high  flow  rates.  This  suggests  that  the  use  of  conven¬ 
tional  vis co-ine r t ial  plot  may  result  in  appreciable 
errors  in  parameter  estimations. 

(5)  The  Klinkenberg  plot  does  not  have  a  unique  shape  but 
is  dependent  upon  the  manner  in  which  the  flow  experi¬ 
ment  has  been  conducted. 

(6)  Since  the  numerical  technique  of  quasilinearization  does 
not  require  data  splitting,  it  is  more  efficient  than 
the  graphical  method. 
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(7)  Values  of  slippage  and  inertial  resistance  coefficients 
obtained  from  published  correlations  serve  adequately  as 
initial  guesses  for  the  quasilinearization  method.  How¬ 
ever,  they  do  not  seem  to  correspond  well  with  the  actual 
values . 

(8)  Confining  pressure  and  the  manner  in  which  it  is  applied 
may  influence  the  results  obtained  in  a  laboratory  flow 
expe  rimen t . 

(9)  The  observation  of  the  same  permeability  values  but  dif¬ 
ferent  slippage  coefficient  and  inertial  resistance  coef¬ 
ficient  values  for  the  same  sample  under  different  experi¬ 
mental  techniques  may  have  been  due  to  shortcomings  in 
either  experimental  or  interpretation  techniques.  How¬ 
ever,  the  possibility  that  this  variation  was  due  to 

some  unexamined  phenomenon  cannot  be  excluded. 

Recommendations 

(1)  Flow  behavior  under  various  constant  downstream  pressures 
should  be  studied  to  investigate  their  effect  on  the  shape 
of  the  Klinkenberg  plot. 

(2)  Steady  state  flow  under  both  varying  upstream  and  downstream 
pressures  should  be  investigated  since  this  more  closely 
simulates  the  conditions  encountered  in  actual  gas  reser¬ 
voirs  . 

(3)  Variations  in  slippage  and  inertial  resistance  coefficients 
under  different  experimental  techniques  should  be  further 
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inves  t igat  ed . 

(4)  The  effect  of  confining  pressure  should  be  studied  in 
greater  detail . 

(5)  If  the  gas  flow  behaviour  at  constant  upstream  pressure 
and  very  low  flow  rates  is  to  be  studied,  more  sophisti¬ 
cated  pressure  regulating  and  measuring  devices  should 

b  e  employed . 

(6)  In  order  to  ensure  uniform  confinement  and  to  avoid 
breakage  of  samples,  the  end  faces  of  samples  should  be 
lapped  as  perfectly  as  possible. 

(7)  Use  of  thick  rubber  gaskets  for  confining  the  cores  should 
be  avoided,  especially  under  higher  pressures.  Otherwise 
the  stress  developed  by  the  flowing  rubber  may  cause  the 
specimen  to  fail. 

(8)  Stress  distribution  in  the  rock  under  confinement  should 
be  examined. 

(9)  The  study  should  be  extended  to  more  competent  rock  samples 
at  pressures  normally  encountered  at  reservoir  conditions. 

(10)  Experimental  verification  of  the  theoretical  expression 
for  inertial  resistance  coefficient  as  given  by  equation 
(13)  should  be  attempted. 

Modification  of  Crafton's  equation  so  as  to  account  for 
molecular  streaming  effects  should  be  attempted. 
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APPENDIX  A 


SALIENT  POINTS  OF  THE  TECHNIQUE  OF 


QUAS IL I NEARI ZAT ION 
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A- 1 


Application  of  the  technique  of  quas i 1 inear i z a t ion 
requires  evaluation  of  the  partial  derivatives  of  the  func¬ 
tional  relationship  with  respect  to  the  parameters. 

Equation  (35)  on  rearrangement  can  be  written  as 
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where  j 


th  e 


t  ll 

j  data  point  and  j  =  1,2,...,N 


t  h 

i  =  the  i  iteration  and  i  =  1,2,„..,M 


f ^  is  the  observed  value  of  f  for  the  data  point 
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From  equation  (A-3) 
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The  nomenclature  used  in  the  computer  program  is 


K 

=  absolute  permeability,  k  (initial  value) 

B 

=  slippage  coefficient,  b  (initial  value) 

FB 

=  inertial  resistance  coefficient,  F  ^  (initial  value) 

THIK 

=  h 

GRAV 

=  G 

EPSLN 

=  £,  the  convergence  criterion 

ROQ 

o 

O' 

it 

XT 

=  T 

PE 

=  pe 

PW 

=  pw 
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etc. 

Additional  nomenclature  used  in  the  computer  output  is 


K1  = 

=  absolute  permeability,  k 

B4  = 

=  slippage  coefficient,  b 

FBI  = 

=  inertial  resistance  coefficient  F 

D 

DCALC(I)  = 

f~^  =  the  function  value  calculated  using  the 

parameters  given  by  quasilinearization. 
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APPENDIX  B 

SAMPLE  CALCULATIONS  USING 


QUASILINEARIZATION  RESULTS 


B  - 1 


Each  of  the  vis  co-inert ial  plots  includes  a  curve 
generated  by  using  the  parameters  estimated  by  means  of 
quasilinearization.  The  procedure  for  generating  such  a 
curve  may  be  illustrated  by  means  of  a  sample  calculation 
as  follows . 

Equation  (35)  which  has  been  derived  to  describe  gas 
flow  through  porous  media  in  radial  systems  may  be  rewritten  as 
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The  application  of  quasi linearization  to  the  data  obtained 
on  sample  111,  in  order  to  estimate  the  parameters  in  the  above 
equation,  resulted  in  the  following: 

k  =  2.719  md 

b  =  12.600  psia 


F 


B 


6.142  x  10 
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ft 


-1 


and 
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Taking  into  account  the  sign  convention  for  ,  equation 
(B-l)  may  be  rearranged  into  the  following  form. 
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If  it  is  assumed  that  the  parameters  obtained  by  means 
of  quasilinearization  are  correct,  they  may  be  substituted  into 
equation  (B-2)  in  order  to  predict  behaviour  under  various 
conditions.  For  example,  given  that 


p^  =  164.52  ps  ia 

Qq  =  0.4432  MS  CF / day 

y  =  0.01759  cp 

r  =  0.48835  ft 
e 

h  =  0.12575  ft 


p  =  13.52  psia 
w 

T  =  526 .511°R 

Z  =0.9985 

r  =  0.08587  ft 
w 


and  that  for  Nitrogen 


G  =  0.9672  (from  NGPSA  Engineering  Data 

Book,  1966,  page  155). 
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one  may  substitute  into  the  right  hand  side  of  equation 
(B-2)  so  that 
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and 
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— - -  =  -  0.08199  (B-5  ) 
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and  the  sum  of  these  three  quantities  is  0.58264. 

Since  this  quantity  is  equal  to  the  left  hand  side  of 

Qo 

equation  (B-2),  it  along  with  the  appropriate  value  of  ( — ) 


may  be  used  to  establish  one  point  of  a  generated  conventional 
vis co-ine r t ial  plot.  The  entire  plot  may  be  generated  by 
repeating  this  procedure  for  other  values  of  the  variables 
involved  » 

Ideally,  this  curve  should  be  identical  to  the  conven¬ 
tional  vis co-inertial  plot  and  deviations  result  from  experi¬ 
mental  error  and  the  fact  that  quasilinearization  fits  the  data 
irrespective  of  their  quality. 
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LISTINGS  OP  COMPUTER  PROGRAMS 
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MAIN 


MAIN 

C  *************************** 


c  * 

C  THIS  PROGRAM  CALCULATES  THE  EXPRESSIONS  NECESSARY  * 
C  TO  MAKE  A  KLINKENBERG  PLOT  IF  THE  VALUE  OF  K  IS  * 
C  READ  IN  AS  EQUAL  TO  1  AND  A  V I  SCO- 1 NER T I AL  PLOT  * 
C  EXPRESSIONS  FOR  K  EQUAL  TO  2  * 
C  IF  K  =  2  9  THE  VALUE  OF  SLIPPAGE  COEFFICIENT  B»  * 
C  SHOULD  BE  READ  IN  TO  GET  THE  EXPRESSIONS  FOR  * 
C  MODIFIED  VISC 0- INERTIAL  PLOT.  IN  CASE  EXPRESSIONS  * 
C  FOR  CONVENTIONAL  V  I  SCO- 1 NERT I AL  PLOT  ARE  DESIRED  * 
C  B  SHOULD  BE  EQUAL  TO  ZERO  * 
C  * 


Q  *************************** 

REAL  KA(50) *  PM ( 50 ) » MW • XX ( 5 ) » ZZ ( 5 ) *  Y  (  5  ) *  X  Y  (  5  ) »PLHS(50) 

*  »QMU ( 50 ) 

READ ( 5  »  7 1  )  N  »  K 

71  FORMAT ( I  3 ♦ 2X  *  1 1 ) 

GO  TO  ( 30  f  3 1 )  iK 

31  R  EAD ( 5  »  3  2  )  B 

32  FORMAT ( F 10 • 5 ) 

30  READ(5t72)RE#RW*THIK 

72  FORMAT ( 3F10 • 5 ) 

WRITE (6 *300) 

300  FORMAT ('  1 «  >20X» • P ( EXT )  P(WELL)  FLOW  RATE  ' 

l'TEMP  VISCOSITY  Z  FACTOR »♦/* 2 IX ♦ 

2 '  ( PS  I A )  (PSIA)  ( MSCF/DAY )  (DEG  F )  ( CP )  •  • / ) 

I  =  1 

301  READ( 5  *  74 ) PE  »PW  *  ROQ  *  T 
74  FORMAT (4F10.5 ) 

*************************** 

* 

EVALUATION  OF  THE  VISCOSITY  AND  COMPRESSIBILITY  AT  * 

THE  ARITHEMATIC  MEAN  TEMPERATURE  AND  PRESSURE  OF  # 

THE  FLOWING  GAS  * 

* 

************************  *  *  * 


XT  =  T+459.7 

AP  =  (PE+PW) /29.3920 

APS=AP^AP 

APC=APS*AP 

XX (  1 ) =504. 

XX ( 2 ) =522. 

XX ( 3 ) =540. 

XX (4) =558. 

XX ( 5 ) =5  76 . 
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APPENDIX  -  C-  2 


MAIN  • • * ( CONT 1 D ) 

MW  =  28.016 

ZZ ( 1 ) =.99998 7E+00-. 3699 8 8E-03*AP+. 20 8731 E-05*APS 
1+.723894E-08*APC 

ZZ ( 2 ) =. 99998 7E+0 0-. 26770 5 E-03*AP+. 2005 99 E-05*APS 
1+.532217E-08*APC 

ZZ ( 3 ) =.99999 5E+00-. 182 03 3E-03*AP+0. 206 590E-Q5*APS 
1+.295211E-08#APC 

ZZ (4) =.999996 E+00-. 102 998 E“03*AP+ . 1 8940 1 E-05*APS 
1+.246034E-08#APC 

ZZ ( 5 ) =.999999 E+00-. 382 71 9E-04*AP+. 1922 74E-05*APS 
1+.705037E-09*APC 

VX= 100.0* ( . 1778E“03*( 1 .+. 8958E-03#( AP-1. ) + . 6 1 2E-0 6* ( AP 

#-l . ) ##2 • 

1+.3997E-07*(AP-1. )#*3. )+.455E-06*( ( ( T-32 . ) /I . 8 ) -2 5 • ) ) 
ZX  =  0.0 

YX= (XT-XX ( 1 ) ) * ( XT-XX ( 2 ) ) * ( XT-XX ( 3 ) ) * ( XT-XX (4 ) ) * ( XT 
*-XX(5)  ) 

DO  410  J=  1  *  5 
W  Y  =  Z  Z ( J ) *YX/ ( XT-XX ( J)  ) 

YDEN= 1 • 

DO  411  M=  1  ♦  5 
IF(J-M)  412*411*412 
412  Y(M)=XX ( J)-XX(M) 

YDEN  =  YDEN*Y ( M  5 
411  CONTINUE 

XY ( J ) =W Y/YDEN 
410  ZX=ZX+XY(J) 

WRITE (6  *302 ) PE  * PW * ROO * T » VX * ZX 

302  FORMAT ( 1 0 1 *18X*6(E10.4*1X) ) 

GO  TO  (600*700) *K 

600  K A (  I  )  =  ( 1424.0*VX*ZX*XT*ROQ*ALOG ( RE/RW ))/((( PE  +  PW ) 

** ( PE-PW ) ) 

1* ( THI K  )  ) 

PM ( I )  =  ( 2*0/ ( PE+PW ) ) 

IF(I-N)  303*304*304 

303  I  =  1+1 
GO  TO  301 

304  WR I TE (6  *77 ) 

77  FORMAT ( ///*5X* 'APPARENT  PERM  ( M. D . ) ' ♦ 5X * ' REC I PROCAl 
*  MEAN  PRESS') 

DO  90  I  =  1  * N 

WRITE (6  *78  )  KA ( I  ) *  PM ( I  ) 

78  FORMAT( 'O' *2( 10X *E13.6) ) 

90  CONTINUE 

700  PLHS(I)  =  ( ( ( ( PE**2 ) - ( PW*#2 ) )  +  ( 2 • *B* ( PE-PW ) ) ) *TH I K ) 
*/  (1424.0 

l*VX*ZX*XT#ROQ ) 

QMU ( I )  =  ROQ/VX 
I F ( I-N )  306  *  306  *  306 
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APPENDIX  -  C 


MAIN 


• . . (CONT’D ) 


305  I  =  1+1 
GO  TO  301 

306  WRITE (6  *80 ) 

80  FORMAT ( / // *  6X  * • LHS  OF  VlSCO  INERTIAL  PLOT«*5X*‘Q 
^/VISCOSITY' ) 

DO  100  I  =  ltN 

WRITE(6  *81 )PLHS(  I ) *QMU( I  ) 

81  FORMAT (  '0»*10X*E13*6*13X*E13.6) 

100  CONTINUE 

800  CALL  EXIT 
END 
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APPENDIX  -  C-  4 

MAIN 

MA  I  N 

C  #**##**#****#■**  #  -«  *#-**#****# 

C  * 

c  #• 

C  TAKING  THE  OBSERVED  FLOW  DATA  AS  INPUT*THIS  PROGRAM  * 

C  GIVES  THE  ESTIMATES  OF  K  *  B  AND  FB  USING  THE  METHOD  * 

C  OF  QUASILINEARIZATION  * 

C  * 

c  * 

C  THE  NOMENCLATURE  AND  UNITS  OF  THE  INPUT  DATA  ARE*  * 

C  RE  =  OUTER  RADIUS  OF  THE  ROCK  SAMPLE  IN  FT .  * 

C  RW  =  BORE  HOLE  RADIUS  OF  THE  ROCK  SAMPLE  IN  FT.  * 

C  THIK  =  THICKNESS  OF  THE  ROCK  SAMPLE  IN  FT.  * 

C  GRAV  =  SPECIFIC  GRAVITY  OF  THE  FLOWING  GAS  * 

C  EPSLN  =  CONVERGENCE  CRITERION  * 

C  K  =  INITIAL  GUESS  FOR  PERMEABILITY  IN  MILL  I D ARC  I ES  # 

C  B  =  INITIAL  GUESS  FOR  SLIPPAGE  COEFFICIENT  IN  PSIA  * 

C  FB  =  INITIAL  GUESS  FOR  INERTIAL  RESISTANCE  * 

C  COEFFICIENT  IN  1/FT.  * 

C  PE  =  PRESSURE  AT  THE  OUTER  RADIUS  IN  PSIA  * 

C  PW  =  PRESSURE  AT  THE  BORE  HOLE  IN  PSIA  * 

C  ROO  =  FLOW  RATE  IN  M$CF/DAY  * 

C  T  =  MEAN  TEMPERATURE  OF  THE  FLOWING  GAS  IN  DEG. F  * 

C  M  =  NUMSER  OF  ITERATIONS  * 

C  N  =  NUMBER  OF  POINTS  IN  THE  DATA  SET  * 

C  COM ( J )  =  APPROPRIATE  COMMENT  WHICH  DESCRIBES  * 

C  THE  DATA  SET  * 

C  * 

C  #  #  *  *  #  #  #  *  #  *  #**f-M-*  # 

c 

£  *****************  #  #  #  #  #  *  #  #  *  * 

c  * 

C  NOTE  -  IN  ORDER  TO  AVOID  CONVERGENCE  TROUBLES  THIS  * 

C  PROGRAM  SHOULD  BE  RUN  IN  DOUBLE  PRECISION  * 

C  * 

Q  #  -ft  *  #  #  -ft  *  #  #  #  #  # 

REAL  D ( 50 )  * P ( 50 )  *R(50)  *S(50) *Q(50) >K*KDET 

REAL  DCALC ( 50 )  *  ERROR ( 50 )  >K1 ( 100 )  *MW 

REAL  B4 ( 100 ) #  FB 1(100)  *  X  X  <  5  )  *  Z  Z  (  5  )  » Y  {  5  )  *XY(5) 

REAL  Cl (50) #C2 ( 50 ) *C3(50) »COM(20) 

READ ( 5  *  1 )  ( COM ( J ) *J  =  1#20) 

1  FORMAT ( 20A4 ) 

WR I TE ( 6  *  2  )  ( COM ( J ) *  J  =  1 *  20 ) 

2  FORMAT (' 0 • *20X * 20A4 ) 

READ ( 5  *  7 1 ) N 

71  FORMAT (13) 
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APPENDIX  -  C-  5 


MAIN  ...(CONT'D) 

C  SETTING  THE  INITIAL  GUESSES 

READ  (5  *500)  K  *  B  *  FB 
500  FORMAT (2F10.5*E13*5) 

READ ( 5*72 ) RE* RW*  TH I K * GR AV ♦ EPSLN 
72  FORMAT ( 5F10. 5 ) 

WR I TE ( 6  * 300 ) 

300  FORMAT (///* 20X ♦’ P ( EXT )  P(WELL)  FLOW  RATE  • 

1 'TEMP  VISCOSITY  Z  FACTOR '  * / ♦  2  IX , 

2 1 ( PS  I A )  (PSIA)  (MSCF/DAY)  ( DEG  F )  (CP)'*/ 

*/) 

I  =  1 

301  READ( 5*74)PE»PW*ROQ*T 
74  FORMAT (4F10. 5 ) 

* 

EVALUATION  OF  THE  VISCOSITY  AND  COMPRESSIBILITY  AT  * 
THE  AR I THEMAT I C  MEAN  TEMPERATURE  AND  PRESSURE  OF  * 
THE  FLOWING  GAS  * 

* 

XT  =  T+459.7 
AP  =  (PE+PW) /29.3920 
APS=AP*AP 
APC=APS#AP 
XX (  1 ) =5  04 • 

XX(2)=522. 

XX ( 3 ) =540 • 

XX (4) =558. 

XX ( 5 ) =576. 

MW=28 .016 

ZZ ( 1 ) =.99998 7E+00-. 36998 8 E-03#AP+. 2087 31 E-05*APS 
1+.723384E-08#APC 

ZZ ( 2 ) =.999987E+00-.267705E-03*AP+*200599E-05*APS 
1+. 53221 7E-08*APC 

ZZ ( 3 ) =.999995E+00-. 182033E-03*AP+0 . 206 590E-05*APS 
1+.295211E-08*APC 

ZZ (4) =.999996E+00- .102 998 E-03*AP+*1 8940 1E-05*APS 
1+.246034E-03*APC 

ZZ ( 5 ) =.999999E+00-.382719E-04#AP+* 1922 74E-05*APS 
1+.705Q37E-09#APC 

VX  =  10C  .0#  (  • 1778E-03*( 1 . + . 8958 E-03* ( AP-1 . )  +  .612E-06* ( AP 
*-l • ) **2 • 

1  +  . 399  7E-0  7* ( AP-1 • >#*3. )+.455E“06#( (  ( T-3  2 .  )  /I • 8 ) -2 5 • )  ) 
ZX  =  0.0 

YX=  (  X  T-XX  (  1  )  )*(XT-XX(2)  )*(XT-XX(3)  )*(XT-XX(4)  )  (  XT 

*-XX ( 5  )  ) 

DO  410  J=  1  *  5 
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APPENDIX  -  C-  6 

MAIN  • • • ( CONT  * D ) 

WY=ZZ ( J)*YX/(XT-XX( J) ) 

YDEN  = 1 • 

DO  411  M= 1  *  5 
IF( J-M)  4 12 >411*412 
412  Y ( M ) =XX ( J)-XX(M) 

YDEN= YDEN#Y ( M ) 

411  CONTINUE 

X  Y  (  J  ) =W Y/ YDEN 
410  ZX=ZX+X Y ( J ) 

WRITE (6  *302 ) PE ♦ P W * ROQ » T ♦ VX  * ZX 

302  FORMAT ( '  '  * 18X *6 ( E10.4 ♦ IX  )  ) 

0(1)  =  ( ( PE**2 ) -PW**2 ) 

Cl ( I )  =  ( 1424.0*VX#ZX*XT*ROQ*ALOG(RE/RW) ) 

C 2 ( I )  =  ( PE-PW ) *2  *0 

C  3 ( I )  =  (0.31602E-11*GRAV*XT*ZX*ROQ*ROQ* ( ( 1 • /RW )  —  ( 1 • 
#/RE) ) ) 

I F (  I-N)  303  *304*304 

303  I  =  1+1 
GO  TO  301 

304  WRITE(6 .77) 

7  7  FORMAT  (///*44X*'K'  *  12X  *  '  B  1  *  12X  »  »  FB  '  ) 

WR I TE ( 6  * 78  ) 

78  FORMAT ( 4 IX  *  '  (MD)  (PSIA)  (1/FT)') 

WRITE(6  *76  )  K  * B  *  FB 

76  FORMAT (' 0 '* 19X *'  INITIAL  GUESS  '  *  4 ( 2 X  * E 1 1 . 4 )  ) 

M  =  1 

26  TPSQ  =  0.0 
TPQ  =  0.0 
TPR  =  0.0 
TOSO  =  0.0 
TOR  =  0.0 
TRSQ  =  0.0 
TPD  =  0.0 
TQD  =  0.0 
TRD  =  0.0 
TPS  =  0.0 
TQS  =  0.0 
TRS  =  0.0 

SETTING  THE  ELEMENTS  OF  THE  MATRIX  AND  THE 
KNOWN  VECTOR 


DO  21  I  =  1 » N 

DCALC(I)  =  (Cl < I )  )/(THIK*K)-(B*C2< I )  )  +  (C3( I )*FB) /(THIK 

#  #  j  (_j  j  i<^  j 

P ( I )  =  -(Cl ( I ) ) / (K*K*THIK) 

0(1)  =  -C2(I) 

R (  I  )  =  (  C  3 (  I  )  )  /  ( THIK*THIK) 

S(I)  =  DCALC ( I ) -K*P (I ) “B*Q ( I ) -FB*R ( I ) 
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APPENDIX  -  C-  7 


MAIN  ...(CONT'D) 

TP  SQ  =  TP  SQ  +P ( I  )  *P( I ) 

TPQ  =  TPQ  +  P ( I )  *  Q ( I ) 

TPR  =  TPR  +  P( I )*R( I ) 

TGSO  =  TQSQ  +  0(1)  *  Q( I ) 

TOR  =  TQR  +  0(1)  *  R ( I ) 

TRSO  =  TRSQ  +  R ( I )  *  R ( I ) 

TPD  *  TPD  +  P ( I )  *  D ( I ) 

TOD  =  TOD  +0(1)  *  D ( I ) 

TRD  =  TRD  +  R ( I )  *  D ( I ) 

TPS  =  TPS  +  P ( I )  *  S ( I ) 

TOS  =  TQS  +  0(1)  *  S ( I ) 

TRS  =  TRS  +  R { I )  *  S ( I ) 

21  CONTINUE 

B 1  =  TP  D-TPS 
B2  =  TQD-TQS 
B3  =  TRD-TRS 

SOLUTION  USING  CRAMER  1 S  RULE 


DDET  =  ( TPSQ# ( TQSQ*TRSQ-TQR*TQR ) ) - ( TPQ* ( TPQ* TRSO -TQR 
**TPR ) ) + ( TPR* ( T 
1PQ*TQR-TQSQ#TPR) ) 

KDET  =  (TPQ*(TOR*B3-TRSQ*B2 ) ) - ( TOSQ* ( T PR*B3“TRS0*B 1 ) ) 
*+( TQR*(TPR*B2 
1-TQR*B1 ) ) 

BDET  =  (TPSQ*(T0R*B3-TRSQ*B2) ) - ( TPQ* ( T PR*B3-TRSG*S 1 ) ) 
*+ ( TPR* ( TPR*B2 
1-T0R*B1 ) ) 

FRDET  =  (TPSQ*(T0SQ*B3-T0R*B2 ) ) - ( TPQ* ( TPQ*B3-TQR*B 1 ) ) 
*+(TPR#(TPQ*B2 
1-T0SQ*B 1 ) ) 

K1(M)  =  KDET /DDE  T 

B4(M)  *  -BDET/DDET 

FBI (M )  =  FBDET/DDET 

WRITE (6  ♦  106 )  M  *  K 1 ( M  > *B4(M) *FB1 (M) 

106  FORMAT ( 20X *' ITERATION  NO . ' *  1 4 ♦ 3 ( 2X ♦ E 1 1 • 4 ) ) 

A  CHECK  ON  THE  CONVERGENCE  OF  THE  SOLUTION 


A 1  =  ABS (  (K1(M)-K)/K) 

A2  =  ABS ( (B4(M)-B)/B) 

A3  =  ABS ( ( FBI ( M ) “FB ) /FB ) 
IF(Al-EPSLN)  82*82*81 

82  IF (A2-EPSLN)  83*83*81 

83  I F ( A3-EPSLN )  6*6*81 
81  K  =  K 1 ( M ) 

B  =  B4 ( M ) 

FB  =  FBI ( M ) 


. 

(nq*( i >q  ♦  * 

(  i )  r  *  (1)0  +  <3:  T  *  C'  T 
) 

(1)2*  (1)9  +  29T  * 

(1)2  *  l ! ) 0  ♦ ■ 2QT  «  20' 

: 

3U/-ITV100  IS 
2  9T • j 9 T  =  19 
2CT-QCT  =  Sa 

„ 2Pt-QPT  «  68 


2  »  r  ■'  *)  ~  •/ .  J 

T ) *V9T ) ♦ ( (99T** 

(  (  •  . ~  ;20T-90T*0fl t 

- 

(  ( I**qQT-I 

S»  #8C  o*?°T)** 
( ( I9*90T-I 

;r«0cT)*paT)+» 

( ( IP*( 20T-X 
fC\T  =  (  )  I 
T -GC\T jCE-  *  ( M I ^ 

;  f  C  )*p,  (»•)  !>•'  (  dOX  *  d )  3T  1 9W 


, 

0 


(>\ ( '- (  ) I> ) ) 28A  *  I A 
-(  ’ )  A 1  ) ) 28A  =  S/ 

9\ ( •  -~ (  )  IS9 ) ) 20A  a  6 a 
{  ,S  *S  (  J21  3-XA)9I 

•  £:  •  <  J2Q3-SAJ3I  SO 

- 

( M) XX  *  x  18 
(  )Aa  =  P 


nnorvnono 


APPENDIX  -  C-  8 


MA  I  N 


. . . (CONT'D ) 


M  =  M  +  l 

IF(M-IOO)  84  *  84  »  6 
84  GO  TO  26 
6  WR I TE ( 6  *  1 08  ) 

108  FORMAT ( »0 •*///// »22X* ' OBS . FUNCT I  ON  CALC. FUNCTION 

*  » 

1 ‘  REL.ERR. ( P • CENT ) « //) 


* 

CALCULATIONS  OF  THE  RELATIVE  PERCENTAGE  ERROR  * 

BETWEEN  THE  OBSERVED  FUNCTION  VALUE  AND  THAT  * 

CALCULATED  USING  THE  PARAMETERS  OBTAINED  # 

FROM  QUASILINEARIZATION  * 

* 


*  *  *  *  *  *  *  ■)(-  #  *  #  #  * 
DO  90  I  =  1 » N 

DCALC(I)  =  (Cl ( I ) ) /(THIK*K )-(B*C2 ( I ) )+ (C3 ( I >#FB) / (THIK 
*#TH I K ) 

ERROR ( I )  =  ABS( ( <D( I )-DCALC< I ) )/D( I ) )*100«0) 

WRITE (6  •  107)  D( I) »DCALC( I ) *  ERROR ( I ) 

107  FORMAT (•  • * 14X # 3 ( 4X »E15 . 5 ) ) 

90  CONTINUE 
CALL  EXIT 
END 


(  •  •  0  D  )  • . . 


* 


,  (  -  )  : 
(fOitd)  ~  !  ; 

( \\  •  (T  1D*^U  •*F3*J3  *1 
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APPENDIX  D 

TABULATED  DATA  AND  RESULTS 


OF  QUASILINEARIZATION 


D-l 


The  following  tables  from  Table  D-l  to  Table  D-8 
contain  the  observed  flow  data.  These  data  are  recorded 
under  a  total  of  9  columns  as  explained  below. 


Column  No . 


Exp 1 ana  t ion 


1 

2 

3 

4 

5 


Serial  number  of  the  observed  data 
point  . 

Upstream  pressure  (p  )  in  inch  of  Hg(*) 
of  psig ( **) . 

Downstream  pressure  (p  )  in  inch  of 
Hg(*)  or  psig  (**)  . 

3 

Time  in  seconds  required  to  flow  50  cm 
(*)  or  1  ft3  (**)  of  gas. 

Inlet  thermocouple  (TC^)  reading  in 
millivol ts . 


6  Outlet  thermocouple  (TC2)  reading  in 

millivolts  . 

7  Barometric  pressure  in  inch  of  Hg. 

8  Room  temperature  in  °F. 


9 


Pressure  shown  by  water  manometer  (WM) 
in  inch  of  water. 


D-2 


TABLE  D- 1 

OBSERVED  DATA  FOR  SAMPLE  111 
Confining  pressure  =  800  psig 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

14.00* 

0.00* 

22 . 84* 

0 .695 

0 .697 

27.71 

67.0 

2 

17.00* 

0,00* 

18.47* 

0 .735 

0 . 741 

• 

CO 

68.5 

3 

14 . 70* 

9  .15* 

53.24* 

0.757 

0 . 759 

27.64 

70.0 

4 

23.50* 

14.89* 

30.82* 

0 . 759 

0 . 762 

27.54 

70.0 

5 

26.99* 

17.20* 

26 . 00* 

0 . 754 

0 . 760 

27.54 

70 .0 

6 

33 . 50** 

0.00* 

1873.50** 

0,757 

0.757 

27.46 

69  .  5 

7 

45.00** 

0.00* 

1206 .02** 

0 . 755 

0 . 755 

27,44 

69  .5 

8 

57.50** 

0.00* 

829 . 77** 

0 .761 

0  .  760 

27.43 

69  .5 

9 

67.00** 

0.00* 

658.00** 

0 . 771 

0.763 

27.43 

69  .5 

10 

80 .50** 

0 .00* 

493 . 50** 

0.770 

0  .  761 

27.43 

69  .5 

11 

92.00** 

0.00* 

398 . 70** 

0 . 755 

0 . 750 

27.53 

69.5 

12 

121.00** 

0.00* 

257 . 22** 

0 . 755 

0  .  751 

27.53 

69 . 5 

13 

151.00** 

0.00* 

176 . 11** 

0.758 

0.751 

27.53 

69  .  5 

14 

182 .00** 

0 .00* 

129 . 14** 

0 . 759 

0 . 751 

27.53 

69  .5 

15 

210 .00** 

0.00* 

101 . 75** 

0  .  761 

0 . 751 

27.53 

69  .5 

16 

239 .00** 

0.00* 

81,56** 

0 . 761 

0  .  751 

27.53 

69  .5 

NB  . 


In  this  run. 


WM  was  not  used 


**00.IS1 

*00.0 

**oo.ees 

1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 
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TABLE  D-2 

OBSERVED  DATA  FOR  SAMPLE  112,  RUN  1 
Confining  pressure  =  500  psig 


2 

3 

4 

5 

6 

7 

8 

9 

4.13* 

0.00* 

15061.42** 

0 .785 

0  .  786 

27.12 

70.0 

0.0 

5 .20* 

0 . 00* 

20.70* 

0 . 754 

0  .  764 

27.36 

70 .0 

0.0 

10.10* 

0 .00* 

10 . 20* 

0 .765 

0 .770 

27.36 

70.0 

0.0 

14 . 80* 

0 .00* 

3716.54** 

0 . 786 

0  c  787 

27.35 

70.0 

0.0 

25.90* 

0.00* 

1914.96** 

0  „  806 

0.807 

27.14 

71.0 

0 . 15 

24.00** 

0 .00* 

843.42** 

0.807 

0 . 806 

27.14 

70  .  5 

0 . 15 

43.00** 

0.06* 

377.95** 

0 . 800 

0 . 800 

27.14 

70.0 

0 . 60 

62.00** 

0.14* 

224 . 41** 

0 . 801 

0 . 801 

27.14 

70.0 

1.00 

72.00** 

0.20* 

179  .  72** 

0 . 790 

0 . 790 

27  .  35 

70.0 

1,30 

84.25** 

0.30* 

141c  73*  * 

0.798 

0 . 798 

27.14 

70.0 

2 . 40 

102 .0** 

0 .53* 

108 , 86** 

0 . 800 

0  o  800 

27 .14 

70 .0 

4 .00 

*00.0 


..  '• 


1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 
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TABLE  D-3 

OBSERVED  DATA  FOR  SAMPLE  112,  RUN  2 
Confining  pressure  =  500  psig 


2 

3 

4 

5 

6 

7 

8 

9 

100 . 0  *  * 

84.0** 

392.13** 

0 . 810 

0 . 806 

27.69 

71.0 

0 .50 

100.0** 

76 .0** 

266  .14** 

0 . 810 

0 . 806 

27.69 

71.0 

0.90 

100 .0** 

67.0** 

206.04** 

0 . 809 

0 . 806 

27.69 

71.0 

1 . 30 

100 .0** 

61 „  0  ** 

184 . 72** 

0.812 

0.80  7 

27 .69 

71.0 

1 . 70 

100 .0** 

52.5** 

161 .42** 

0.810 

0 . 805 

27.69 

71.0 

2 . 10 

100.0** 

42.5** 

144 . 36** 

0 . 810 

0 . 807 

27.69 

71.0 

2 .50 

100 .0** 

35.0** 

135 . 70** 

0 . 810 

0.807 

27.69 

71.0 

2  .  80 

100 .0** 

26.0** 

126.97** 

0 . 809 

0 . 806 

27.69 

71.0 

3 . 10 

100 .0** 

18.0** 

122 .05** 

0 . 810 

0 . 810 

27.69 

71.0 

3.  70 

100 .0** 

12.0** 

119 .09** 

0.812 

0 . 812 

27.69 

71.0 

3.90 

100 .0** 

0 .53* 

113.78** 

0.812 

0 . 813 

27.69 

71.0 

4.10 

x 

**0.001 
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TABLE  D-4 

OBSERVED  DATA  FOR  SAMPLE  112,  RUN  3 
Confining  pressure  =  500  psig 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

80.0** 

75 .4** 

1427 . 82** 

0 .777 

0  .  779 

2  7.60 

71.5 

0.05 

2 

80 .0** 

58.5** 

351 . 71** 

0 . 786 

0.785 

27.60 

71.5 

0.50 

3 

80.0** 

46.0** 

246  .  72** 

0 . 796 

0  .  792 

27.60 

71.5 

1.00 

4 

80.0** 

34.5** 

204  .'72** 

0.802 

0  .  798 

27.60 

71.5 

1 . 30 

5 

80.0** 

23.0** 

178.48** 

0.807 

0 . 801 

27.57 

71 . 5 

1.60 

6 

80 .0** 

13.0** 

165 . 83** 

0 . 809 

0 . 800 

27.57 

71.5 

1.90 

7 

80.0** 

0 . 10** 

156.69** 

0 . 814 

0 . 804 

27.57 

71.5 

2 . 10 

I  i  ‘ 


7* 

' 

1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 
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TABLE  D-5 

OBSERVED  DATA  FOR  SAMPLE  2,  RUN  1 
Confining  pressure  =  500  psig 


2 

3 

4 

5 

6 

7 

8 

9 

200.0** 

0.0* 

1154 . 86** 

0 . 800 

0  .  796 

27.67 

68.00 

0 . 15 

180 .0** 

0.0* 

1422 . 36** 

0 . 804 

0 . 799 

27.67 

68.00 

0 . 15 

160 .0** 

0.0* 

1792.13** 

0 . 806 

0 . 801 

2  7.67 

68.00 

0 . 10 

140.0** 

0.0* 

2326 .51** 

0.808 

0.803 

27.67 

68.25 

0.05 

120 .0** 

0.0* 

3120 .21** 

0 . 805 

0 . 805 

27.67 

68.50 

0.03 

100 .0** 

0.0* 

4409 . 45** 

0 . 805 

0 . 805 

27.65 

68.50 

0 .0 

80.0** 

0.0* 

6582.68** 

0 . 795 

0  .  795 

27.65 

69 .00 

0 . 0 

60.0** 

0.0* 

10582 .68** 

0 . 804 

0 . 800 

27.65 

69.00 

0.0 

4  8.0*  * 

0.0* 

26.40* 

0 . 800 

0  .  798 

27.66 

6  8.6 

0 . 0 

36.0** 

0.0* 

40 . 50* 

0 . 796 

0 . 795 

27.66 

6  8.6 

0 .0 

28.0** 

0.0* 

58 . 35* 

0.793 

0.791 

27.66 

6  8.6 

0 . 0 

22.0** 

0.0* 

82  .2* 

0 . 795 

0  .  794 

2  7.66 

68.6 

0.0 

16.0** 

0.0* 

122.4* 

0 . 790 

0 . 790 

27.64 

70.0 

0.0 

8 .0** 

0 .0* 

286.95* 

0 . 790 

0  .  790 

27.66 

6  8.6 

0.0 

t 

.  • 
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TABLE  D-6 

OBSERVED  DATA  FOR  SAMPLE  2,  RUN  2 
Confining  pressure  =  500  psig 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

200.0** 

183.0** 

12.60* 

0 .833 

0 . 825 

27.69 

70.0 

0.0 

2 

200.0** 

160.0** 

3425 .20** 

0 . 840 

0 . 830 

27 . 69 

70.0 

0.0 

3 

200.0** 

142.0** 

2543.31** 

0 . 835 

0 . 830 

27.69 

70.0 

0.08 

4 

200.0** 

108.0** 

1815 .22** 

0 . 800 

0 . 802 

27.69 

69  .2 

0.10 

5 

200.0** 

80.6** 

1530 .71** 

0 . 795 

0 . 800 

27.69 

69 .0 

0.10 

6 

200 .0** 

56 .0** 

1385.83** 

0.792 

0.798 

27.69 

68.5 

0 .15 

7 

200.0** 

30.0** 

1289  .76** 

0 . 790 

0 . 795 

27.69 

68.5 

0 . 15 

8 

200.0** 

14.0** 

1244.09** 

0 . 790 

0.795 

27 . 69 

68 . 3 

0.15 

9 

200.0** 

0 .0* 

1217.85** 

0 . 790 

0.794 

27.69 

68 . 3 

0 .15 

i 


■ 

*'  £  ■  ,  c  :  x 

**0 . 00£ 

**90.  <>■  :>£I 

1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 
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TABLE  D-7 

OBSERVED  DATA  FOR  SAMPLE  2,  RUN  3 
Net  Confining  pressure  =  400  psig 


2 

3 

4 

5 

6 

7 

8 

9 

10 .0** 

0.0* 

204.90* 

0 . 845 

0.847 

27.71 

71.50 

0.0 

16.0** 

0 .0* 

114.60* 

0 . 854 

0 . 852 

27.73 

71 . 60 

0.0 

20.0** 

0.0* 

84.90* 

0 . 834 

0 . 831 

27.73 

70.50 

0.0 

22.0** 

0.0* 

75.00* 

0 . 840 

0 .841 

27.71 

71.30 

0.0 

28.0** 

0.0* 

54 .60* 

0 . 835 

0 . 831 

27.74 

70.50 

0.0 

40.0** 

0.0* 

32.40* 

0 .836 

0 . 840 

27.71 

71.30 

0.0 

60.0** 

0.0* 

17.1* 

0 . 830 

0 . 831 

27.71 

71.00 

0  .0 

80.0** 

0.0* 

6137 . 54** 

0 . 812 

0 . 812 

27.80 

69  .  70 

0 .02 

100.0** 

0.0* 

4230.37** 

0.813 

0 . 815 

27.80 

69  .  70 

0.03 

120.0** 

0.0* 

3105.38** 

0 . 817 

0 . 817 

27 . 80 

69  .  70 

0.05 

140.0** 

0.0* 

2393. 70** 

0 . 817 

0 . 819 

27 . 80 

69  .  70 

0.05 

160.0** 

0.0* 

1891 . 86** 

0 . 823 

0 . 821 

27 . 80 

69 .90 

0.07 

180 .0** 

0.0* 

1541.19** 

0 .823 

0 . 823 

27 . 80 

70.00 

0.10 

200 .0** 

0.0* 

1281 . 89  ** 

0 . 825 

0.824 

27 . 80 

70 .00 

0.10 

**?8  ,I8£I 

**0. 00£ 

/ 
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TABLE  D- 8 

OBSERVED  DATA  FOR  SAMPLE  2,  RUN  4 
Net  Confining  pressure  =  400  psig 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

200.0** 

0 .0* 

1259  .  84** 

0 . 787 

0 .785 

27.48 

68.0 

0 .15 

2 

200 .0** 

14 . 38** 

1286 .17** 

0.796 

0  .  791 

27.48 

68.0 

0 . 15 

3 

200 .0** 

28.63** 

1314.96** 

0 . 796 

0 . 794 

27.48 

68.0 

0 . 10 

4 

200 .0** 

52 . 88** 

1401.58** 

0 . 831 

0 . 821 

27.64 

69.6 

0.0 

5 

200 .0** 

87  .  75** 

1623.62** 

0 .836 

0.826 

27.64 

69.6 

0.0 

6 

200.0** 

120 .0** 

2047.24** 

0 . 841 

0 . 831 

27.64 

69.5 

0.0 

7 

200 .0** 

133.50** 

2363  .  70** 

0 . 873 

0 . 876 

27.68 

73.2 

0 .0 

8 

200 .0** 

165.95** 

4183.83** 

0 .880 

0 . 882 

27.65 

74.0 

0.0 

D 
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TABLE  D  -  9 

PROCESSED  DATA  AND  QUAS I L I NE AR I Z AT  I  ON  RESULTS  SAMPLE  111 

P(EXT)  P(WELL)  FLOW  RATE  TEMP  VISCOSITY  Z  FACTOR 


(PSIA) 

(PSIA) 

(MSCF/DAY) 

(DEG  F) 

(  CP  ) 

0 

•  2048  E 

02 

0 

•  1 3  6  1 E 

02 

0 

•6100E-02 

0 

•  6389E 

02 

0 

•1745E-01 

0 

•  9996E 

00 

0 

•2203E 

02 

0 

•1368E 

02 

0 

•7560E-Q2 

0 

• 6577E 

02 

0 

. 1749E-01 

0 

•9996E 

00 

0 

• 2079E 

02 

0 

•  1 8  06E 

02 

0 

•  2  590E-02 

0 

•  6668  E 

02 

0 

• 1752E-01 

0 

•  9996E 

00 

0 

•2506E 

02 

0 

•2Q83E 

02 

0 

•4460E-02 

0 

• 6678E 

02 

0 

•1753E-01 

0 

•9996E 

00 

0 

•  2  678  E 

02 

0 

•  2 1 97E 

02 

0 

•5290E-02 

0 

•  6663  E 

02 

0 

•1752E-Q1 

0 

•  9  99  5  E 

00 

0 

•  469  8  E 

02 

0 

• 1348E 

02 

0 

•4284E-01 

0 

•  666 1 E 

02 

0 

•1753E-01 

0 

•  9  9  9  4  E 

00 

0 

•5847E 

02 

0 

•  1 347E 

02 

0 

•6450E-01 

0 

•  66  53  E 

02 

0 

•  1753  E-0 1 

0 

•  9993E 

00 

0 

•7097E 

02 

0 

•  1 347E 

02 

0 

•9372E-01 

0 

• 6679E 

02 

0 

•1755E-01 

0 

• 9993E 

00 

0 

• 8047E 

02 

0 

•  1 347E 

02 

0 

•  1 18  IE  00 

0 

•6709E 

02 

0 

•1756E-Q1 

0 

•  9992E 

00 

0 

•9397E 

02 

0 

•  1 347E 

02 

0 

•  1 57 5E  00 

0 

•  6702  E 

02 

0 

•1757E-01 

0 

•  999 1 E 

00 

0 

•  10  5  5  E 

03 

0 

• 1352E 

02 

0 

•1957H  00 

0 

•  6643  E 

02 

0 

•  1 756E-C 1 

0 

•  9990E 

00 

0 

•  1 345  E 

03 

0 

•  1 3  5  2  E 

02 

0 

• 3034E  00 

0 

• 6644E 

02 

0 

•1757E-01 

0 

•  9987E 

00 

0 

• 1645E 

03 

0 

•  1 3  5  2E 

02 

0 

•  443  IE  00 

0 

•  66  5 1 E 

02 

0 

•1759E-01 

0 

•9985E 

00 

0 

•1955E 

03 

0 

• 1352E 

02 

0 

•  6 04 3 E  00 

0 

• 6653E 

02 

0 

•1761E-01 

0 

•  9  9  8  3  E 

00 

0 

•  2235E 

03 

0 

• 1352E 

02 

c 

•7670E  00 

0 

• 6659E 

02 

0 

•1763E-01 

0 

•  998  1 E 

00 

0 

•  2  52  5  E 

03 

0 

•  1 3  5  2E 

02 

0 

•9569E  00 

0 

•  66  57E 

02 

0 

•1764E-01 

0 

• 9979E 

00 

K 

B 

FB 

(MD) 

(PSIA) 

(1/FT) 

INITIAL  GUESS 

0.4679E 

01 

0.6255E 

01 

0.1679E 

09 

ITERATION 

NO* 

1 

0.1307E 

01 

0.1259E 

02 

0.6145E 

10 

ITERATION 

NO. 

2 

0. 1986E 

01 

0. 1259E 

02 

0.6144E 

10 

ITERATION 

NO. 

3 

0.2522E 

01 

0. 1259E 

02 

0.6145E 

10 

ITERATION 

NO- 

4 

0.2705E 

01 

0.1259E 

02 

0 • 6 1 45  E 

10 

ITERATION 

NO. 

5 

0.2719E 

01 

0. 1259E 

02 

0.6145E 

10 

ITERATION 

NO. 

6 

0.2719E 

01 

0.1259E 

02 

0 • 6 145E 

10 
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TABLE  D  -  9 


OBS. FUNCTION  CALC • FUNCT I  ON  REL • ERR .( P • CENT ) 


0 

• 23445E 

03 

0 

•  2  3022E 

03 

0 

•18050E 

01 

0 

•  2  98  3  1 E 

03 

0 

*  2  9  29 1 E 

03 

0 

.  18069E 

01 

0 

•10594E 

03 

0 

•  10420E 

03 

0 

•  1 644 1 E 

01 

0 

•19414E 

03 

0 

.  19136E 

03 

0 

•  1 43 1 2  E 

01 

0 

•  2  3444E 

03 

0 

•23206E 

03 

0 

•  10169E 

01 

0 

•  20258E 

04 

0 

•20266E 

04 

0 

•37512E-01 

0 

• 32379E 

04 

0 

•  3  1966E 

04 

0 

•12746E 

01 

0 

•48555E 

04 

0 

• 48672E 

04 

0 

• 24076E 

00 

0 

•  62943E 

04 

0 

•63033E 

04 

0 

•14273E 

00 

0 

•  8  6493E 

04 

0 

•  86654E 

04 

0 

•  1 8692  E 

00 

0 

•  1 095 1 E 

05 

0 

. 10988E 

05 

0 

• 33058E 

00 

0 

•17913E 

05 

0 

•  17787E 

05 

0 

•  70353E 

00 

0 
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05 
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05 

0 

•  5  2448  E 

00 
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•38045E 

05 
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05 
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0 

•49779E 

05 
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05 
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00 
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TABLE  D  -  10 

PROCESSED  DATA  AND  QUAS I L  I  NE AR I Z AT  I  ON  RESULTS  SAMPLE  112  RUN  1 


P(EXT) 

P(WELL) 

FLOW  RATE 

TEMP 

VISCOSITY  Z 

FACTOR 

(PSIA  ) 

(PSIA) 

( MSCF/DAY ) 

(DEG  F) 

(CP) 

0.1534E 

02 

0.1332E 

02  0.5100E-02 

0.6793E 

02  0.1755E-01 

0.9997E 

00 

0.1599F 

02 

0.1343E 

02  0.6610E-02 

0.6653E 

02  0.1751E-01 

0.9997E 

00 

0.1839E 

02 

0.1343E 

02  0.1355E-01 

0.6711E 

02  0.1753E-01 

0.9997E 

00 

0.2070E 

02 

0.1343E 

02  0.2090E-C1 

0 • 6797E 

02  0.1755E-01 

0.9997E 

00 

0.2605E 

02 

0.1332E 

02  0.4010E-01 

0 .6888E 

02  C.1758E-01 

0.9996E 

00 

0.3732E 

02 

0.1332E 

02  0.9113E-01 

0 .6888E 

02  C.1758E-01 

0.9995E 

00 

0.5632E 

02 

0.1335E 

02  0 • 2037E  00 

0 .6859E 

02  0.1758 E -01 

0.9994E 

00 

0.7532E 

02 

0.1339E 

02  0*3 440E  00 

0 .6863E 

02  0 • 1 760E-0 1 

0.9992E 

00 

0.8343E 

02 

0.1353E 

02  0.4327E  00 

0.6813E 

02  0# 1759E-01 

0.9992E 

00 

0#9757E 

02 

0.1347E 

02  0.5461E  00 

0.6850E 

02  0.1760E-01 

0  #9991E 

00 

0.1153E 

03 

0. 1359E 

02  0.7141E  00 

0  •  68 59E 

02  0.1762E-01 

0.9989E 

00 

K 

B 

FB 

(MD) 

(PSIA)  (1/FT) 

INITIAL 

.  GUESS 

0.1661E  02 

0  •  38 1 5  E 

01  0 • 9235E 

09 

ITERATION 

NO#  1 

0.1081E  02 

0 • 4562  E 

01  0 • 5  547E 

10 

ITERATION 

NO.  2 

0.1213E  02 

0 • 4562  E 

01  0.5547E 

10 

ITERATION 

NO#  3 

0# 123  IE  02 

0 • 4562  E 

01  0.5547E 

10 

ITERATION 

NO#  4 

0 • 12 3  IE  02 

0 • 4562  E 

01  0.5547E 

10 

ITERATION 

NO#  3 

0.1231E  02 

0 • 4562  E 

01  0.5547E 

10 

•  • 


•  CONTD 


nufl  SIX 


0 -  C  3J8AT 

>T JJ2  ?  MOITASIfAM.  IJieAUC  QUA  AT AQ  O-22'Or^ 


t*)a3  %  nuooaiv 

(Q-)  (3  3i0) 


3TAP  .‘CJ3 
( YA0\"O?*  ) 


(  jj'wjq 
( a : i 


( TX3 ) 9 
( AI2^) 


3T 9 •  C  I C-3c  c  Y I • -  >0 

:  Y  9  ?  9  •  C  1 0-312’'  I  •  - 

•  X  i  *3 £  '  • 

.  -'-jScYI* 

.  -  • 

•  ,  10—3  0  :?\  •  ■> 

,0  i  -  i  8  $  Y  £  •  S  '■■) 

. 

. 

<v  9.  I  -  S  d  Y  I  •  0  SO 


3??Yd*  SC-3001 S*C  S; 

,  s  ■"  -  *  0  X  s  a  •  .  S 

.  ~ 

.  -  • 

•  i  •-  g  ;  •' • 

3  xo-“i€ r $ 

-3  p  ^  h  a  #  0  3  Y  S  •  0  S 

!€d3d.O^  30^C#0  S< 
£lPd*G  (0  3YS£A# 

c  e  p  a  •  ‘  g o  .ii  •**<■  • 


S££X#3  SO  3A£cI»C 
^£[#0  SO  3W«tO 
3e A£  X  #0  SC  I • 

3£^£l #0  SO  30Y0S* 
3S££I*0  SO  3eOdS.O 
1 S £ £  l • 0  SO  3 S £  Y £  •  0 
1 £  £  £  I • 0  SC  3  S  £  d  c  •  " 
19EfcX#0  SO  3S££r«0 
•  e  e  £  1  •  S'  3EAe°#G 
lYAti.C  SO  3YPYP.C 
39e£X»0  £0  3 £ e 1 1 • 


(  T^Vi  > 


9  j  e  t' s  y  . 

3 Y  Ac8  • 

-  f  j  Y  c  C  •  G 

cl  3Y**Se» 


i a; >q) 


(0  ) 


ic  3Sa5A# 

io  osae^# 
s  a  t.u  • 
sacA. 


s  Ur  ox# : 

: 

s  .  oexs:  .c 

su  3 x £ s x • o 

SO  3l£SI#0 

uesx.o 


saauo  JAiTini 
x  #on  HO  I  T  A913T  I 
.0/  nCITAc^I 
.on  noiTAflBTi 
^  .On  no  I  T  AflBT  i 
.On  '0ITAP3TI 


•  • 


D  -  13 


TABLE  D  -  10 


OBS • FUNCT  I  ON  C ALC • FUNCT I  ON  REL • ERR .( P . CENT  ) 


0.58153E 

02 

0.75164E 

02 

0.15793E 

03 

0.24813E 

03 

0.50096E 

03 

0.12158E 

04 

0 • 29945E 

04 

0  •  54950E 

04 

0.71157E 

04 

0 • 93402E 

04 

0.13116E 

05 

0 • 60048  E 

02 

0.78095E 

02 

0. 16356E 

03 

0.25762E 

03 

0*51196E 

03 

0.12359E 

04 

0.29937E 

04 

0.54390E 

04 

0.71089E 

04 

0  •  94097E 

04 

0. 130-90E 

05 

0  •  32574E 

01 

0.38999E 

01 

0.35685E 

01 

0.38264E 

01 

0  •  2 1969E 

01 

0.16538E 

01 

0 • 2  8469E-0 1 

0*101 89E 

01 

0.94417E-01 

0.74453E 

00 

0.19511E 

00 

01  -  C  -3JAAT 


0  • ! 
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389005.0 

38eeax.  ; 

40 

a^eesi *o 

40 

385ISX.O 

I0« 

—  >848 s. 

40 

3Yt'P9S.O 

40 

3c  499£ • 0 

398I0X.Q 

4lh309£45.0 

40 

305945. 0 

-3YI449.0 

40 

39  80  X  V  #0 

40 

1V5XXV.0 

3f544Y«C 

40 

3Y9Q49.0 

40 

3SC459.0 

3X  I  e°x#c 

eo 

30908X  .0 

eo 

38XXOX.O 
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TABLE  D  -  11 

PROCESSED  DATA  AND  QUAS I L I  NEAR  I Z AT  I  ON  RESULTS  SAMPLE  112  RUN  2 


0 


P(EXT) 

P ( WELL ) 

FLOW  RATE 

TEMP 

VISCOSITY  Z 

FACTOR 

(PSIA  ) 

(PSIA) 

(MSCF/DAY) 

(DEG  F) 

(CP) 

•1136E 

03 

0.9760E 

02 

0.2203E 

00 

0.6895E 

02 

0.1767E-01 

0.9984E 

00 

•  1 1 36  E 

03 

0  •  8  9  6  0  E 

02 

0.2949E 

00 

0 .6895E 

02 

0.1767E-01 

0.9984E 

00 

•  1 1 36  E 

03 

0.8060E 

02 

0.3814E 

00 

0.6893E 

02 

0.1766E-01 

0.9985E 

00 

•  1 1 36  E 

03 

0.7460E 

02 

0.4258E 

00 

0 .6902E 

02 

0. 1766E-01 

0.9985E 

00 

•  1 136E 

03 

0 • 66 1 OE 

02 

0.4878E 

00 

0.6893E 

02 

0.1765  E-0 1 

0.9986E 

00 

•1136E 

03 

0 • 56 1 OE 

02 

0  •  546  IE 

00 

0 .6897E 

02 

0.1 765E-0 1 

0.9987E 

00 

•  1 1 36  E 

03 

0.4860E 

02 

0.5814E 

00 

0.6897E 

02 

0.1764E-01 

0.9987E 

00 

•  1 1 36  E 

03 

0 .3960E 

02 

0.6218E 

CO 

0 .6893E 

02 

0.1764E-01 

0.9988E 

00 

•1136E 

03 

0.3160E 

02 

0.6479E 

00 

0 .6904E 

02 

0. 1764E-01 

0.9988E 

00 

•1136E 

03 

0.2560E 

02 

0.6643E 

00 

0-6913E 

02 

0.1 764E-0 1 

0.9989E 

00 

•  1 1 36  E 

03 

0 • 1 3  8  6E 

02 

0.6971E 

00 

0.6915E 

02 

0. 1763F-01 

0.9990E 

00 

K 

B 

FB 

(MD) 

( PSIA ) 

(1/FT) 

INITIAL 

GUESS 

0. 1661E 

02 

0.3815E 

01 

0  •  9235E 

09 

ITERATION 

NO.  1 

0. 1065E 

02 

0.1398E 

02 

0.9151E 

10 

ITERATION 

NO.  2 

0. 1202E 

02 

0.1398E 

02 

0.9151E 

10 

ITERATION 

NO.  3 

0. 1222E 

02 

0. 1398E 

02 

0.9151E 

10 

ITERATION 

NO.  4 

0. 1223E 

02 

0.1398E 

02 

0.9151E 

10 

ITERATION 

NO.  5 

0. 1223E 

02 

0.1398E 

02 

0.9151E 

10 
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TABLE  D  -  11 


OBS*  FUNCT ION 

CALC. FUNCT 

ION 

REL.ERR. ( 

P.C 

0 

•  3  3  79  2  E 

04 

0.34606E 

04 

0.24097E 

01 

0 

•48768E 

04 

0  «  47728  E 

04 

0.21322E 

01 

0 

•64086E 

04 

0  •  64  340  E 

04 

0.39732E 

00 

0 

•  7  3  39  8  E 

04 

0.73083E 

04 

0.42924E 

00 

0 

•  85357E 

04 

0.85822E 

04 

0 • 5  44  8  8  E 

00 

0 

•  97  577E 

04 

0  •  97929E 

04 

0.36086E 

00 

0 

* 10543E 

05 

0.10517E 

05 

0.23768E 

00 

0 

•11336E 

05 

0. 11364E 

05 

0.24303E 

00 

0 

•  1 1906E 

05 

0. 11871E 

05 

0.29647E 

00 

0 

•  1 2249E 

05 

0.12170E 

05 

0.64721E 

00 

0 

•  1  2  7 1 2  E 

05 

0. 12781E 

05 

0. 54066E 

00 
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TABLE  D  -  12 

PROCESSED  DATA  AND  QUAS I L I NE AR I Z AT  I  ON  RESULTS  SAMPLE  112  RUN  3 


P(EXT) 

P  ( WELL ) 

FLOW  RATE 

TEMP 

VISCOSITY  Z  FACTOR 

( PS  I  A  ) 

(PSIA) 

( MSCF/DAY ) 

(DEG  F) 

(CP) 

0  •  93  5  5E 

02 

0.8895E 

02 

0  •  5  4  6  2  E  ■ 

-01 

0 

•  67  59E 

02 

0.1762E-01 

0.9985E 

00 

0.9355E 

02 

0.7205E 

02 

0 • 2  2  2  OE 

00 

0 

•  6793E 

02 

0.1762E-01 

0.9987E 

00 

0.9355E 

02 

0.5955E 

02 

0.3169E 

00 

0 

•  6  8  3 1 E 

02 

0 • 1 762  E-0 1 

0.9988E 

00 

0.9355E 

02 

0.4805E 

02 

0.3822E 

00 

0 

•6859E 

02 

0.1762E-01 

0.9988E 

00 

0.9354E 

02 

0.3654E 

02 

0.4383E 

00 

0 

•  68  77E 

02 

C • 1762E-01 

0 • 9989E 

00 

0.9354E 

02 

0.2654E 

02 

0.4721E 

00 

0 

•  68  79E 

02 

0. 1762E-01 

0.9990E 

00 

0.9354E 

02 

0.1364E 

02 

0.4999E 

00 

0 

• 6900E 

02 

Q.1762E-01 

0.9991E 

00 

K  '  B  FB 


(MD) 

(PSIA) 

(1/FT) 

INITIAL  GUESS 

0.1661E 

02 

0.3815E 

01 

0 • 92  3  5  E 

09 

ITERATION 

NO. 

1 

0. 1076E 

02 

0.2097E 

01 

0 . 4344E 

10 

ITERATION 

NO. 

2 

0. 1210E 

02 

0.2096E 

01 

0 • 4344E 

10 

ITERATION 

NO. 

3 

0.1228E 

02 

0  *  2097  E 

01 

0.4344E 

10 

ITERATION 

NO. 

4 

0.1229E 

02 

0.2097E 

01 

0 • 4343  E 

10 

ITERATION 

NO. 

5 

0. 1229E 

02 

0.2098E 

01 

0.4344E 

10 

ITERATION 

NO. 

6 

0. 1229E 

02 

0.2096E 

01 

0.4344E 

10 

•  • 
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TABLE  D  -  12 


OBS. FUNCTION 

CALC • FUNCT 

ION 

REL • ERR • (P 

•  CE 

0  #  83955E 

03 

0.83805E 

03 

0.17824E 

00 

0.35606E 

04 

0#  35680E 

04 

0.20658E 

00 

0.52058E 

04 

0 • 52223E 

04 

0.31816E 

00 

0 .64433E 

04 

0 • 63986E 

04 

0  •  6932 5 E 

00 

0.74146E 

04 

0 • 74329E 

04 

0.24607E 

00 

0.80455E 

04 

0.80550E 

04 

0.11816E 

00 

0.85638E 

04 

0.85589E 

04 

0.57853E- 

'01 

: . 

ipeaos.c 
oo  3 aieie.o 
)0  3eS£?a.O 

. 

cj  3ei8i;. 
XO- J£29T£.C 


3?0e£8«0,.  £0  3229£8#0 

3  0  0  d  2  £  •  ^  AO  3d0d2£.0 

;£ssse.o  *o  33eos2.o 

A-  ~:d89£3»0  AO  3££AAd •  0 

A  39S-CAT.-C  AO  38AI AT •  0 

BCceCt'.O  AO  322A09.0 

AO  398228t0  AO  33£d2B.O 
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TABLE  D  -  13 

PROCESSED  DATA  AND  QUAS I L I  NEAR  I Z AT  I  ON  RESULTS  SAMPLE  2  RUN  1 


P(EXT) 

P(WELL) 

FLOW  RATE 

TEMP 

VISCOSITY  Z 

FACTOR 

(PSIA  ) 

(PSIA) 

( MSCF/DAY ) 

(DEG  F) 

(CP) 

0 

•  2  1 35  E 

03 

0.1359E 

02 

0 

•6817E-01 

0.6850E 

02 

0.1 767E-0 1 

0.9982E 

00 

0 

•  1935E 

03 

0.1359E 

02 

0 

. 5  535E-0 1 

0 . 6  8  6  5  E 

02 

0.1766E-01 

0.9984E 

00 

0 

•  1 73  5  E 

03 

0.1359E 

02 

0 

•4392E-01 

C .6875E 

02 

C.1765E-01 

C.9985E 

00 

0 

•  1 5  3  5  E 

03 

0.1359E 

02 

0 

•3381E-01 

0 .6884E 

02 

C.1764E-01 

0.9987E 

00 

0 

•  1335E 

03 

0  •  1359E 

02 

0 

•  2  5  20E-0 1 

0.6381E 

02 

0.1763E-01 

0.9988E 

00 

0 

•  U35E 

03 

0.1358E 

02 

0 

•  1 78QE-0 1 

0 .6881E 

02 

0.1762E-01 

C.9990E 

00 

0 

•  93  58  E 

02 

0.1358E 

02 

0 

•1192E-01 

0 • 6836E 

02 

0.1 760E-0 1 

0 • 999  IE 

00 

0 

•7358E 

02 

0.1358E 

02 

0 

.7  380E-02 

0  *  68 68 E 

02 

0.1 760E-0 1 

0.9993E 

00 

0 

•  6  1 5  8  E 

02 

0.1358E 

02 

0 

•  5  2  50E-02 

0.6854E 

02 

0.17  5  9E-0 1 

C.9993E 

00 

0 

•  49  5  8  E 

02 

0.1358E 

02 

0 

. 3420E-02 

0.6838E 

02 

0.175  8  E-Q 1 

0.9994E 

00 

0 

•4158E 

02 

0.1358E 

02 

0 

•2370E-02 

0 .6822E 

02 

0.1 75  7E-0 1 

0.9995E 

oc 

0 

•  3558E 

02 

0.1358E 

02 

0 

•  1 680E-02 

0 .6834E 

02 

0.1 757E-0 1 

0.9995E 

00 

0 

•2957E 

02 

0.1357E 

02 

0 

.1 120E-02 

0 .6813E 

02 

0.1756E-01 

0.9996E 

00 

0 

•  2  1 5  8  E 

02 

0  •  1  3  5  8  E 

02 

0 

•  4  800E-03 

0.6813E 

02 

0. 1755E-01 

C.9997E 

00 

K 

B 

FB 

(MD) 

(PSIA) 

( 1/FT  ) 

INITIAL 

GUESS 

0. 

89 76 E  00 

0. 1043E 

01 

0.3429E 

11 

ITERATION 

NO.  1 

0. 

6746 E  00 

0.9196E 

00 

-0.7662E 

11 

ITERATION 

NO.  2 

0. 

7162E  00 

0.9196E 

00 

-0.7661E 

11 

ITERATION 

NO.  3 

0. 

7139E  00 

0.9196E 

00 

-0 . 7662  E 

11 

ITERATION 

NO.  4 

0. 

7 1  89 E  00 

0.9196E 

00 

-0.7661E 

11 
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TABLE  D  -  13 


OBS. FUNCTION  CALC . FUNCT I  ON  REL . ERR •( P . CENT ) 


0 

•  4  543  5  E 

05 

0 

•45329E 

05 

0 

•  23400E 

00 

0 

•  37292E 

05 

0 

#  37400E 

05 

0 

•28967E 

00 

0 

•  2  9948E 

05 

0 

•  30083E 

05 

0 

•  448  2  5  E 

00 

0 

•  2  340  5  E 

05 

0 

•  2  3424E 

05 

0 

•  8  39  37E-0 1 

0 

•  1 766 1 E 

05 

0 

.  17612E 

05 

0 

•27861E 

00 

0 

•12716E 

05 

0 

•  1 2  5  2  2  E 

05 

0 

•  1 5 1 97E 

01 

0 

•  8  572  8  E 

04 

0 

• 84069E 

04 

0 

•  1 93  57E 

01 

0 

• 52296E 

04 

0 

•52184E 

04 

0 

•  2  1338E 

00 

0 

•  3608  IE 

04 

0 

• 37096E 

04 

0 

•28116E 

01 

0 

•  2  2  74 1 E 

04 

0 

•24113E 

04 

0 

•60350E 

01 

0 

•15447E 

04 

0 

. 16664E 

04 

0 

•78732E 

01 

0 

•  10817E 

04 

0 

•11786E 

04 

0 

•  895  36E 

01 

0 

•  6904 1 E 

03 

0 

•  78  3 17E 

03 

0 

.13434E 

02 

0 

•28136E 

03 

0 

•  3  3  37  5  E 

03 

0 

•  1 86 1 8  E 

02 

91  -  Q 


El  -  JJPA- 


OtOJ.5  a«J3/ 

WO  IT 

:  ’U3.  -j/ o 

W0ITDWU3.280 

300A£S#0 

3ae*a4.o 

ri  v  d  e  b  S  •  0 

30i  4TC.0 

3S9ST&.0 

3as?A+**  :» 

"SO Of .0 

1  -499S.0 

10- 

•3TE9€  .0 

'  A  A  ■  .  C 

3204£S.0 

3ia=i^^.c 

ISidYX  .0 

eo 

3XddVX«0 

10 

arcidi.o 

3$sasx.o 

aXTSX." 

-vavv'i.o 

3  9  a  0  ^  6  •  0 

4Q 

ips^ee.o 

3 v  £  ?  X  S  •  •' 

4  x 

48xse.o 

40 

3d9SS?.0 

10 

Boll  £*C 

40 

3*90T€.0 

40 

3 1 80d£  •  0 

3exi4£*o 

.40 

X4TSS.0 

sevsv.o 

i4dddX  *0 

40 

3  V 445  X  •  0 

10 

3dCa9B«0 

idPTXX.O 

3TX0OX*O 

3-  i**€X*0' 

3VX€8V.C 

X4C90#O 

3*  I  331 .0 

;aV£  £.0 

3d£X8S.O 
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TABLE  D  -  14 

PROCESSED  DATA  AND  QUAS I L  I  NEAR  I Z AT  I  ON  RESULTS  SAMPLE  2  RUN  2 


P(EXT) 

P ( WELL ) 

FLOW  RATE 

TEMP 

VISCOSITY  Z  FACTOR 

(PSIA) 

(PSIA) 

( MSCF/DAY ) 

(DEG  F) 

(  CP  ) 

0.2136E 

03 

0  •  1966E 

03 

0*1 099E-0 1 

0.6987E 

02 

0. 1780E-01 

0.9971E 

00 

0.2136E 

03 

0.1736E 

03 

0.2291E-01 

0.7013E 

02 

0.1780E-01 

0.9973E 

00 

0.2136E 

03 

0.1556E 

03 

0.3085E-01 

0.7002E 

02 

0 • 1 779E-0 1 

0 • 9974E 

00 

0 • 2 1 36E 

03 

0.1216E 

03 

0.4330E-01 

0.6863E 

02 

0.1773E-01 

0.9975E 

00 

0.2136E 

03 

0 • 942  OE 

02 

0.5136E-01 

0 .6847E 

02 

0*1771 E-0 1 

0.9977E 

00 

0.2136E 

03 

0.6960E 

02 

0.5680E-01 

0  •  6836E 

02 

0.1770E-01 

0.9978E 

00 

0.2136E 

03 

0.4360E 

02 

0*6 10  3E-0 1 

0.6825E 

02 

0. 1768E-01 

0.9980E 

00 

0.2136E 

03 

0.2760E 

02 

0.6329E-01 

0.6825E 

02 

0.1767E-01 

0.9981E 

00 

0 • 2 1 36E 

03 

0.1360E 

02 

0 • 6466E—0 1 

0.6822E 

02 

0.1766E-01 

0.9982E 

00 

K 

B 

FB 

(MD) 

( PSIA ) 

(1/FT) 

INITIAL  GUESS 

0.8976E 

00 

0. 1043E 

01 

0.3429E 

11 

ITERATION 

NO* 

1 

0.6857E 

00 

0.3595E 

02 

0.2717E 

12 

ITERATION 

NO* 

2 

0.7239E 

00 

0.3595E 

02 

0.2717E 

12 

ITERATION 

NO* 

3 

0  •  72  6 1 E 

00 

0*  3595E 

02 

0.2717E 

12 

ITERATION 

NO* 

4 

0*  726  IE 

00 

0.3595E 

02 

0.2717E 

12 

ITERATION 

NO. 

5 

0.7261E 

00 

0.3595E 

02 

0.2717E 

12 
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TABLE  D  -  14 


OBS. FUNCTION  C ALC • FUNCT I  ON  REL . ERR .( P • CENT ) 


0.69734E 

04 

0.15488E 

05 

0.21413E 

05 

0.30838E 

05 

0.36731E 

05 

0  •  40780E 

05 

0 • 43724E 

05 

0.44863E 

05 

0.45439E 

05 

0.71793E 

04 

0*  15501E 

05 

0 • 2 1 326  E 

05 

0.30674E 

05 

0 • 36895E 

05 

0.40870E 

05 

0.43578E 

05 

0.44925E 

05 

0*  45440E 

05 

C.29535E 

01 

0.86504E* 

-01 

0.40828E 

00 

0.53148E 

00 

0.39244E 

00 

Q.21966E 

00 

0 • 332  88  E 

00 

0.13964E 

00 

0.82634E- 

-03 

*1  -  0  1JFAT 


(  T'/3D.C1 5  •  *  *J:* 


0 1 TD/U3 • 0 J  ‘O 


lO!TDHin»2*0 
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00  3SS8QA.0 
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00  3 3 d 9 1 S  *  - 
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00  j*d?£X*0 


3E9VIV.0 
Hioeei .c 
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C  3*VdG€#Q 
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4^  )  3^£Y93«0 

ao  3B?^ei«o 

eO  3£IAIS*0 
60  38EB0E  •  0 

eo  3ieTd£.o 

•!C  308Y0A.0 

ao  3AST£A«0 
ao  3£d°A^«0 
ao  3<?£AaA«o 
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TABLE  D  -  15 

PROCESSED  DATA  AND  OUAS I L I  NEAR  I Z AT  I  ON  RESULTS  SAMPLE  2  RUN  3 

P(EXT)  P ( WELL )  FLOW  RATE  TEMP  VISCOSITY  2  FACTOR 


(PSIA) 

(PSIA) 

(MSCF/DAY) 

(DEG  F) 

(CP  ) 

0 

•  2  360E 

02 

0.1360E 

02 

0 

•6700E-03 

0 

•  7063E 

02 

0 

•1762E-01 

0 

•  9997E 

00 

0 

•  2  96  1  E 

02 

0.1361E 

02 

0 

•1200E-02 

0 

•  7095  E 

02 

0 

• 1763E-01 

0 

•  9996E 

00 

0 

•  3362E 

02 

0.1362E 

02 

0 

•  1 630E-02 

0 

•  7002  E 

02 

0 

•1761E-01 

0 

•  9  996  E 

00 

0 

•  3  560E 

02 

0.1360E 

02 

0 

•  1 840E-02 

0 

• 7038E 

02 

0 

•1762E-01 

0 

•  9996E 

00 

0 

•  4 16  2  E 

02 

0  •  1 3  6  2  E 

02 

0 

•  2  530E-02 

0 

•7004E 

02 

0 

•1761E-01 

0 

•  9995E 

00 

0 

• 5360E 

02 

0.1360E 

02 

0 

.4260E-02 

0 

• 7027E 

02 

0 

•1763E-01 

0 

• 9994E 

00 

0 

• 7360E 

02 

0. 1360E 

02 

0 

•8090E-02 

0 

•  6993  E 

02 

0 

•1763E-C1 

0 

•  9993E 

00 

0 

•9365E 

02 

0.1365E 

02 

0 

•1284E-01 

0 

•  69  1 3  E 

02 

0 

•  1 76  2  E-0 1 

0 

•  9  99 1 E 

00 

0 

•  1 1 36E 

03 

0.1365E 

02 

0 

•1663E-01 

0 

•  6921 E 

02 

0 

•1763E-01 

0 

• 9990E 

00 

0 

• 1336E 

03 

0.1365E 

02 

0 

.2538E-01 

0 

• 6934E 

02 

0 

•1765E-01 

0 

•  9  9  8  8  E 

00 

0 

•1536E 

03 

0.1365E 

02 

0 

•3293E-01 

0 

• 6939E 

02 

0 

•1766E-01 

0 

•  9987E 

00 

0 

• 1736E 

03 

0.1365E 

02 

0 

•4 165E-01 

0 

•  69  56  E 

02 

0 

.1767E-01 

0 

•  9  9  8  5  E 

00 

0 

• 1936E 

03 

0.1365E 

02 

0 

•5112E-01 

0 

• 6960E 

02 

0 

.1769E-01 

0 

•  9984E 

00 

0 

•  2 136E 

03 

0.1365F 

02 

0 

•6146E-01 

0 

• 6967E 

02 

0 

•1770E-01 

0 

•  998  3E 

00 

K 

B 

FB 

(MD) 

(PSIA) 

(1/FT) 

INITIAL  GUESS 

0.6976E 

00 

0. 1043E 

01 

0.3429E 

11 

ITERATION 

NO. 

1 

0.6470E 

00 

0.6884E 

01 

0 • 46  5  2  E 

11 

ITERATION 

MO. 

2 

0.6974E 

00 

0.6884E 

01 

0.4653E 

11 

ITERATION 

NO. 

3 

0.7016E 

00 

0.6884E 

01 

0.4653E 

11 

ITERATION 

NO. 

4 

0.7017E 

00 

0.6884E 

01 

0.4652E 

11 
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TABLE  D  -  13 


OBS. FUNCTION  CALC . FUNCT I  ON  REL • ERR .( P . CENT ) 


0 

•  372  1 9E 

03 

0 

•  36543E 

03 

0 

•  1 8 1 72  E 

01 

0 

•69133E 

03 

0 

•68223E 

03 

0 

•  1  3863E 

01 

0 

•94498E 

03 

0 

•  947  27  E 

03 

0 

•24184E 

00 

0 

•10528E 

04 

0 

. 10791E 

04 

0 

•  33657E 

00 

0 

•  1 5469E 

04 

0 

•  1 5 139  E 

04 

0 

•  2 1370E 

01 

0 

•26887E 

04 

0 

•  2  6  548  E 

04 

0 

•12637E 

01 

0 

•  5  2  3  3 1 E 

04 

0 

•  5  2  74  1 E 

04 

0 

•78189E 

00 

0 

•  85846E 

04 

0 

•85918E 

04 

0 

•  83765E 

-01 

0 

•  1 273  OE 

05 

0 

.  12756E 

05 

0 

•  1 9890E 

CO 

0 

•  1 7676E 

05 

0 

•  17710E 

05 

0 

•19184E 

00 

0 

•23423E 

05 

0 

•  2  3  344E 

05 

0 

•  33556E 

00 

0 

•  2  9969E 

05 

0 

•  29989E 

05 

0 

•68909E 

-01 

0 

•37315E 

05 

0 

•  37  3 1 3  E 

05 

0 

•  87655E 

-02 

0 

•  4546 1 E 

05 

0 

•45467E 

05 

0 

•  1 22  2  5  E 
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TABLE  D  -  16 

PROCESSED  DATA  AND  QUAS I L I  NEAR  I Z AT  I  ON  RESULTS  SAMPLE  2  RUN  4 


P(EXT) 

P(WELL) 

FLOW  RATE 

TEMP 

VISCOSITY  Z  FACTOR 

(PSIA) 

(PSIA) 

( MSCF/DAY ) 

(DEG  F) 

(CP) 

0  •  2  1 34E 

03 

0.1349E 

02 

0 

•6206E-C1 

0 -6795E 

02 

0.1765E-01 

C  •  9  9  8  2  E 

00 

0.2134E 

03 

0.2737E 

02 

0 

•6079E-01 

0.6829E 

02 

0 • 1 767E-0 1 

0  •  998 1 E 

00 

0.2134E 

03 

0.4212E 

02 

0 

•5945E-01 

0 .6836E 

02 

0. 1768E-01 

0.9980E 

00 

0.2135E 

03 

0.6645E 

02 

0 

•5592E-01 

0.6973E 

02 

Q.1773E-01 

0 • 9979E 

00 

0.2135E 

03 

0 • 1 0 1 3E 

03 

0 

•4827E-01 

0.6995E 

02 

0*1775  E-Q 1 

0.9977E 

00 

0.2135E 

03 

0.1335E 

03 

0 

•  3  8  29E-0 1 

0 .7018E 

02 

0*1778  E-Q 1 

0.9975E 

00 

0  •  2  1 3  5  E 

03 

0.1470E 

03 

0 

•3298E-Q1 

0.7189E 

02 

0* 1783E-01 

0.9975E 

00 

0  •  2  1 3  5  E 

03 

0 • 1 79  5E 

03 

0 

•1858E-01 

0.7218E 

02 

0. 1785E-01 

0.9974E 

00 

K 

B 

FB 

(MD) 

(PSIA) 

( 1/FT ) 

INITIAL  GUESS 

0.8976E 

00 

0.1043E 

01 

0.3429E 

11 

ITERATION 

NO. 

1 

0.7231E 

00 

0.3061E 

01 

0.7835E 

11 

ITERATION 

NO. 

2 

0.7504E 

00 

0.3060E 

01 

0.7835E 

11 

ITERATION 

NO. 

3 

0.75 14E 

00 

0.3065E 

01 

0.7835E 

11 

ITERATION 

NO. 

4 

0.7514E 

00 

0.3062E 

01 

0  *  7836E 

11 

ITERATION 

NO. 

5 

0 • 75  1 5  E 

00 

0.3059E 

01 

0.7836E 

11 

ITERATION 

NO. 

6 

0.7515E 

00 

0.3063E 

01 

0.7835E 

11 

ITERATION 

NO. 

7 

0.7515E 

00 

0  *  3062  E 

01 

0.7836E 

11 

iteration 

NO. 

8 

0.7514E 

00 

0.3062E 

01 

0.7836E 

11 
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TABLE  D  -  16 


OBS*  FUNCT ION  C ALC  .  FUNCT I  ON  REL 


0.45398E 

05 

0 • 4 5  548  E 

05 

0. 

0 • 44804E 

05 

0.44674E 

05 

0. 

0.43806E 

05 

0.43704E 

05 

0. 

0.41198E 

05 

0.41228E 

05 

0. 

0.35347E 

05 

0.35399E 

05 

0« 

0.27772E 

05 

0.27832E 

05 

0. 

0.23985E 

05 

0.23973E 

05 

0. 

0.13385E 

05 

0. 13308E 

05 

0* 

•ERR. (P.CENT ) 


32923 E  00 
28906E  00 
23387E  00 
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49930E-01 
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